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A UNIVERSAL PNEUMATIC MULTIPLIER=DIVIDER AND A 
SQUARE~ROOT EXTRACTOR 


Iu. I. Ivlichev and E, M. Nadzhafov 


(Moscow) 


Block schematics are given for a universal pneumatic multiplier-divider and a 
square-root extractor. The equations describing the operation of these devices are 
derived and estimates of the accuracy of these devices are given. The results of 
experimental investigations are also cited, 


In automating production processes, the necessity arises for implementing the operations of multiplication, 
division and square~root extraction on signals. The use of pneumatic devices for carrying out these operations is 
advantageous since, as will be made clear in the sequel, such devices are both simple and reliable, 


The present paper is devoted to the description and analysis of a universal pneumatic multiplier-divider 
and a square-root extractor, both developed in the Pneumo-hydro-automation Laboratory of the Institute for 
Automation and Remote Control of the AN SSSR. 


The multiplier-divider permits execution of the operation: 


P,P, 
= (1) 


where P, is the output pressure, and Py, P, and Pgs are input pressures, All pressures are expressed in “excess” units 
(excess dver atmospheric pressure). By holding one or two of the inputs fixed, it is possible to execute the follow- 
ing additional simple operations: 


1) multiplication of two signals: Pg = kPyP3; 

2) squaring of a signal: Py = kP*; 

3) division of one signal by another: P, = kP;/ Pg; 

4) division of the square of one signal by another signal: Py = P*/ Ps; 
5) amplifying a signal within a wide range of gain variation: Py = kPy; 
6) obtaining the reciprocal of a quantity: Py = k/P; 


In all the foregoing formulas, k is a scale factor which may be varied within broad limits. All the pressures 
(input and output) are kept within the range, 0 to 1 atmosphere. Therefore, any quantity within this range can 
be multiplied or squared, Division, multiplication by a constant factor k, and obtaining the reciprocal must be 
limited, due to the constraint that the output magnitude not exceed one atmosphere. 


Equation (1) shows that the device is comprised of an equalizing bridge circuit which provides the relation- 
ship: Pg/P, = P,/Ps. In executing the computational functions implemented by the multiplier-divider, the 
condition P, <P, will hold. For the implementation of the universal operation given in (1), the condition, 

Ps = Py miust be opserved, 


If necessary, one can broaden the domain of the mathematical functions implemented by a preliminary 
division of each of the variables Py, P, and Ps, by a previously established quantity k in an ordinary pneumatic 
divider, in the form of a flow chamber with valves at the input and the output. One of these valves is controlled, 
allowing the magnitude of k to be varied. 


If, in (1), the output quantity Py is set equal to Ps then, formally, one can obtain the oniee of square- 
root extraction, expressed by the equation: 


Ps = 2, (2) 


where, in this case, Pg is the output quantity of the square-root extracter. In practice, obtaining the square root 
in this way is difficult and, therefore, another circuit, based on the scheme for the multiplier-divider, is ‘ 
for the realization of (2). 


If one of the two input quantities, Py or P,, is fixed, it is then possible to extract the square root of one sig- 
nal with a necessary scale factor introduced, i.e., the operation Ps = k/P, where k = /Py or k = /Pg. The scale 
factor, k, can be varied within wide limits, If the operation in (2) is implemented, the condition Ps< P, will 
hold, 


1. Description of the Multiplication and Division Block 


If, in a nozzle-valve type relay (Fig. 1), there is a simultaneous variation of the input pressure, P;, (before 
the constant valve) and the position of the flapper of the variable valve, then the pressure, P,, in the inter-valve 
chamber will depend on these two variable quantities, In the general case, P, is a nonlinear function of some 
one of these two quantities, the other being fixed, even if a linear dependency is maintained between the discharges 
from the variable and constant valves and the pressure drops across them. If the nozzle opening is controlled by 
the other external pressure, P,, then under certain conditions (we shall call them “compensation conditions"), it 
is possible to realize the function: 


Px = kPsP,. (3) 


This relationship underlies the universal pneumatic multiplier- 
divider we have developed. 


EDP The schematic of this device is shown in Fig. 2. The device 
consists of two pneumatic relays, the nozzles of which are overlapped 
by one valve flapper. Input pressures P, and Ps are supplied to the 
constant valves of these relays, and input pressure P; is applied to the 
ciosed chamber blocked by a membrane. The blocking membrane 
functions like a servo system, In chamber B there is always developed 
a pressure, P,, equal to the pressure, Py, in chamber A, With a variation 
in pressure Py, the position of the piston of the blocking membrane is 
changed, executing the role of the valve flapper, and the pressure, 
Py, of the inter-valve chamber, B, changes until equilibrium with the 
: other valve opening is attained, This is quaranteed by the inflexible 
4 R, a, negative feedback which transmits the pressure P, from the inter-valve 
: chamber of the relay to chamber B of the blocking membrane. If 
4 pressure P, is constant, but pressure Ps across the constant relay valve 
Fig. 2. changes, the cited pressure equilibrium of chambers A and B is attained 
by a change in the nozzle opening. The equilibrium states, deter- 
mined by the equality of pressures Py and P,, are reached by small movements of the piston (about 0.05 mm). 
Thus, either pressure P, or pressure Ps controls the nozzle which acts as a valve on the output of air from chamber 
B. Since the butt end of the piston controls the opening of the two symmetrically placed nozzles, pressures Py 
and Py also control the nozzle aperture of output chamber C. The pressure in the inter-valve chamber C of the 
second relay thus depends both on the nozzle aperture (and, consequently, on pressures Py and Pg) and on pressure 
P, at the input of the constant valve. The output pressure of the multiplier-divider is pressure Py. An increase 
of pressures P, and P, engenders an increase in pressure Py, while an increase in pressure Ps decreases pressure P,. 
If the compensation conditions are realized for the two pneumatic relays, controlled by a common valve flapper, 
then a functional dependency of the form given in (1) can be obtained. 
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2. Static Properties of the Multiplier-Divider 


In order to develop the static properties of this block, we set up the equations for the steady-state processes 
in the flow chambers of the pneumatic relays and the equation for equilibrium of the forces on the piston. 


We shall assume as given a linear dependency of the mass discharge on the pressure drops across all the 
valves. If this condition holds, then the steady-state processes in the pneumatic relays can be described by the 
equations: 


Ry (4) 


(5) 


where R is the pneumatic impedance of a valve, defined as R = AP/G. Here, AP is the pressure drop across the 
valve and G is the mass discharge per second from the valve. 


If we ignore small parameters (weights of the moving parts, jet reaction and membrane elasticity) and 
assume equality of the effective areas of the membrane, we get, from the condition of equilibrium of forces on 
the piston, 

Py =P, (6) 


By taking into account the fact that the section of the variable valve's effective aperture is proportional 
to the flapper’s travel, we can determine the variable valve’s impedance by means of an expression which is in 
good correspondence with experimental results: 


h h 
Ra = R' =, Rea = R" =", (7) 


where h is the flapper's travel measured from the nozzle, hy, is the maximum working travelof the flapper,and 
R® and R® are the inherent impedances of the nozzles. 


The system of equations (4) through (7) describes the static properties of the multiplier-divider. Solving 
this system with respect to the output quantity Py gives 


— 1) (8) 


To obtain the necessary functional relationship (1), it is necessary that the following compensation con- 
dition hold, 


R 
(9) 


which implies the requirement that the valve impedances of the pneumatic relays be proportional. In particular, 
this condition holds if the valve impedances are equal: Ry = Rg, R’ =R°. 


The realization of condition (9) in practice is hindered by the circumstance that, in the case when the 
relay impedances are proportional, the nozzles would have to function on proportional, and not equal, gaps, This 
creates a difficult design problem. We shall therefore consider the compensation condition (9) in the particular 
case of equal impedances, 


Thus, the device whose schematic is shown in Fig. 2.realizes the operation Py = PyP,/Py if, in the steady~ 
state, there is a linear dependence of the discharges from all valves, constant and variable, on the pressure drops 
across them, and if compensation condition (9) holds. The linearity of the discharges and the compensation con- 
dition are attained by the proper choice of impedances and by providing the nozzles with wide butt ends. 
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3. Static Error of the Multiplier-Divider 


If the conditions given above held exactly, then the operation Py = PyP,/ Pg would be executed without error. 
Actually, errors arise for three reasons: 


a) due to deviations of the discharge characteristics of the valves from the linear; 
b) due to the inexact meeting of compensation condition (9); 


c) due to the effect of the small factors which were ignored in setting up Eq. (8) (weight of the moving 
parts, jet reaction, membrane elasticity, etc). 


We now determine the effects of these causes on the accuracy of the device's functioning, 


a) Estimate of the error arising from deviations of the discharge characteristics from linearity. In writing 
Eq. (4) and (5), we assumed that, for all valves, the discharge depended linearly on the pressure drop, Actually, 
this is not the case, and for two reasons: the nonlinearity of the impedance and the compressibility of air. To 
take these factors into account, we introduce nonlinear terms into the equations for the steady=state processes 
in the pneumatic relays: 


P, 
Ri + $1 (Ps, Px) = Ra + (Px), (10) 
+ $2 (Ps, Pa) = (Ps), (11) 


where %3. %g» cq ANd gz are terms which take into account the deviations of the corresponding valve discharges 
from the linear. By taking Eq. (6) and (7) into account, we get 


R. 
P 
where a = R° Ry/R'Rg. 
To determine the relative error, we put expression (12) in the form: 
Rie (1 +8,), (13) 


(a —1)(1— $1) eq) — — 
1+ (@—1) +p 


We now use compensation condition (9), and, for simplicity, we add the assumption that Ry = Ry =R. Asa 
result, we get 


(14) 


Ps (2 — cq) — Pa ($1 — 
(FF +91 — 


(15) 


It is clear from this that a small relative error can be obtained if the conditions of nonlinear compensation 
hold, 


= Pca Fi = Pei» (16) 


which fs difficult to achieve in practice. To obtain an estimate of the error, we assume that the differences of 
the deviatiors from linearity of the discharges are equal for both relays, and have the maximum value: 


(1 — Pc1) = (P2 — Fea) = Ag. (17) 


Then 


P; de) 
By introducing the notation Ps/R =C, we obtain a formula for the relative error in the form 


Ae 


14° 


For sufficiently small Ag/G, one can consider that 


= 


5, < . (19) 
For example, for capillaries of 0.22 mm diameter and 40 mm length, nozzle diameters of 0.5 mm and 
nozzle butt end diameter of 6 mm, this error does not exceed 0.2%, 


b) Estimate of the error due to the inexact meeting of compensatién condition (9). We assume that com- 
pensation condition (9) does not hold exactly, and that there is an initial error of 


i" 


1. 
With this, Eq. (7) will have the form: 


P _ (1+ PiP2 
+ 


The relative error due to the inexactitude of the compensation is determined from the expression: 


With a small initial error, €,, the greatest relative error, 5g, obeys the relationship 
(21) 


In practice, the error, 53, can be made as small as desired, 
c) Estimate of the error connected with the neglected small parameters. If the small parameters (piston 
mass and jet reaction) are taken into account, Eq. (6) for the equilibrium of forces on the piston takes the form: 


where Fer is the total effective membrane area, Q is the mass of the moving parts attached to the piston and 
Ncy, Ncg are the jet reaction forces of the nozzles, The elasticity of the membrane was ignored in setting up 
Eq. (22), since the working travel of the valve is very small. 


If the effective nozzle areas are denoted by F,4 and Fog, the jet reactions can be estimated by the in- 
equalities: 


Na<PsFa, Nea< P2F ca 
We rewite Eq. (22) for the case when the jet reactions assume their maximum values; 


| 
| 
| 
PaFage + Q— Nex — Nea = 0, (22) i 
PyFegq— + Q— — = 0. (23) q 
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If we take Eqs. (4), (5) and (7) into account, we get 


(24) 


where a has the same denotation as in Eq. (12). For the initial force error we take the expression 


Feff 
If compensation condition (9) holds (a = 1), Eq. (24) takes the form: 


(25) 


Pym (26) 


The relative error related to the small parameters obeys the condition: 


€2 
(27) 
With P, = 1 in the case when the weight of the moving parts is 50 grams, the effective total membrane area 
is 20 cm* and the nozzles have the dimensions given above, this error will be 0.15%. 
The total static relative error can be defined as the sum: 


bz = 5, + 8,48, (28) 
and is-about 0.35%, If one of the pressures serves as a fixed scale factor, another error may arise due to the 
instability of this pressure. 

4. Certain Questions of Multiplier-Divider Dynamics 


For estimating certain dynamic properties of the device, we will derive its equation, under the assumption 
that the condition of static compensation (9) holds; the weight of the moving parts attached to the piston can 
be ignored due to its smallness, and the motion of the piston occurs in the working region of the nozzle gap. * 


With these assumptions, the processes in the flow chambers as described by the equations: 


(29) 
Py Pe GPs 30 


where Cy = V/RT is the total capacity of the flow chamber and the feedback chamber, C, = V,/RT is the total 
capacity of the output flow chamber and the load, V, and Vz are the volumes of the corresponding chambers, R 
is the gas constant and T is the absolute temperature. 


If Eqs. (6) and (7) are joined to Eqs, (29) and (30), then the differential equation of the device's dynamics 
can be obtained: 


PT, + P,[a + —a) = (31) 


where T; = RyCy and T, = R,C, are the time constants of the flow chambers. With a = 1, the equation takes the 
form: 


d 
+ Py (Ps = P,P,. 


*In its open-looped state, the pneumatic relay has a gain on the order of 4000, and is shunted by an inflexible 
negative feedback path. 


Q— PsF., — Pol, 
|| 
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As is obvious from Eq. (31), the functional relationship of the input and output pressures in the device, even 
in the simplest model, is described by a complicated differential equation. Even if the static compensation con- 
dition holds and the equation is thereby simplified, it still remains sufficiently complex, Equation (32) goes over 


into a linear differential equation with constant coefficients only in the particular case of device operation with 
Py = const and P, = const, 


5. Description of the Universal Device for Square-Root Extraction. 


As was mentioned above, P, must be set equal to Ps in Eq. (1) to obtain the operation of square-root extraction. 
However, due to the circumstance that the pneumatic relay, shunted by a feedback path, loads the second relay, 
the necessary accuracy can not be obtained. In this case it would be necessary to satisfy the compensation con- 
dition (9), not on the basis of an equality of the valves’ 


iui | impedance, but on the basis of a choice of differing 
a> valves satisfying the compensation condition and 
wr | LOI decreasing the influence of the load. Such a solution 
a7 | would be very complex and difficult for practical 
; Re, nin Ree = realization. It is much simpler to execute the operation 
of square-root extraction by introducing some decoupling. 
The block schematic for the square-root extractor 
Rh is shown in Fig. 3. The device operates the same as the 
1 multiplier, with the condition that Pg = Py. This equality 


indicates the presence of a negative feedback path from 

Fig. 3. the device's output to the input of the first relay. The 
presence of this feedback has the effect that the output 
quantity itself delays an increase of its value. 


If the compensation condition holds, the circuit given for extracting square roots executes operation (2) 
with a definite accuracy. 


6. Static Error of the Square-Root Extracter. 


If the conditions assumed in the analysis of the multiplier-divider were to hold exactly then, if there were 
no internal load, the operation Ps = PP, would be executed without errors. 


In the square-root extractor, errors arise from the following causes; 

a) from the influence of the inherent load; 

b) from the deviation of the valves’ discharge characteristics from the linear; 
c) from the inexact meeting of compensation condition (9); 


d) from the effect of small factors which were neglected in setting up Eq. (8) (weight of the moving parts, 
jet reaction, membrane elasticity, etc). 


We now estimate the effect of these factors on the accuracy of the square-root extractor’s operation. 


a) Estimate of the error due to the effect of the inherent load. The condition of equilibrium of the forces 
acting on the membrane servo system is defined by the equality of the air pressure above and below the mem- 
brane: 


Ps = Py: (33) 


If there is just one output valve, the opening of which is controlled by the membrane, and if this valve is 
capable of providing a sufficiently small pressure to the submembrane chamber, then for any input valve there 
is attained an ideal addition of the inter- valve chamber pressure to the controlling pressure. It is necessary for 
this that the supply pressure be not less than the controlling pressure, 


If there is anadditional uncontrolled valve discharging air from the intervalvye chamber, then this additional 
condition is not always fulfilled. 


q 
aa . 
4 
ij 
a 
al 
\ 
‘a 
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| 
fi 
fi 
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For a small value of controlling pressure, the additional discharge from the intervalve chamber does not 
entail an error end the condition of pressure equivalence is observed,due to the variation of the discharge nozzle 
in such a manner that a balance of intake and discharge in the intervalve chamber occurs, corresponding to 
equivalence of the values of the pressures. 


For the limiting value of controlling pressure, at which the discharging nozzle is completely closed, the 
air pressure in the intervalve chamber will be less than the controlling pressure. This difference is the larger, the 
smaller the impedance of the second output valve. 


In the servo system scheme considered, the presence of the feedback path imposes an additional condition, 
the result of which is that the greatest difference between the controlling pressure and the intervalve chamber 
pressure occurs for a value of controlling pressure less than the limiting one. However, it is not difficult to see 
that the error, in the case considered, is always less than the greatest error which occurs with an open nozzle. 


To estimate the error, we consider an isolated servo system with a maximal load. In this case, the equation 
for the steady-state process in the servo system's flow chamber will have the form: 


Ps (34) 


where Rg is the impedance of the supply valve, Res is the impedance of the servo system's nozzle and Py is the 
supply pressure. 


The exact value of the output pressure, Ps , corresponding to the absence of a load, is determined from 
Eq. (34) by setting Ry = ©. With this, we get 


(35) 


a Po 36 
P, = (36) 
1+ Res hh 


If 44 is a small quantity, it then follows from Eq. (37) that 


Ry Ro 
| A, (1 + (38) 


In choosing the valve impedances for the flow chamber of the servo system, it is necessary to add to con- 
dition (38) an additional condition, related to the necessity of eyacuating the servo system's flow chamber at 
maximum travel of the valve flapper. This condition is obtained from the equality of the discharges from the 
flow chamber when the nozzle is open and when the given residual chamber pressure is Ps, For simplicity we 
ignore the load, thus obtaining 

Pyo—P, P, 


= 39 
i (39) 


where R™ is the inherent impedance of the open nozzle. 
Denoting the relative magnitude of the residual pressure by 0, = P§/P 9, we get 


1 
1 + 
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1+ 
If there is a load, it then follows from Eq. (34) that 
and the relative error is determined from the estimate: 
R 
41 - (37) 
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If oy is small, we may write 


Rr” 


The greatest error occurs for Reg = R™. 
By substituting (41) in (38), we get 
R 
Rs | 4; (1 ++ (42) 


Thus, by taking conditions (41) and (42) into account, we 
can, given one of the impedance ratios, find the other if the 
admissible error A, is given. 


| 


In pneumatic computing devices, an effective choice of 
the magnitude of oy is of the order of 0,01. 


Without given valve dimensions, the error related to the 
presence of loading does not exceed 0,25% The remaining 
errors occurring in the square-root extractor are about half the 
corresponding errors in the multiplier-divider. 


The total static relative error can be determined as the 
sum of all the errors, and amounts to 0.5%, 


1. Equation of the Square-Root Extractor 
Dynamics. 


We set up the equation for the square-root extractor 
dynamics under the assumptions that all static errors equal zero 
and that the piston motion occurs within the working range of 
the valve gap. The remarks made in the first paragraph of 
section 4 should also be kept in mind, With these assumptions, 
the dynamic equation of the square-root extractor may easily 
be obtained from dynamic equation of the multiplier-divider, 
(32), by adjoining to it the equations describing the dynamic 
processes in the servo systems. 


With no inherent load, the processes in the servo system 


are described by the equation 
Po—Ps dP, 43 


where Cs = Vs/RT is the total capacity of the load and the 
following chamber, Vs is the volume of the corresponding 
chambers, Res is the impedance of the discharging nozzle: 


Rg = (44) 


where z is the working travel of the servo system's membrane 
and z,, is its maximum working travel. 


The equilibrium equation for the servo system's mem 
brane. has the form: 


(Py — Ps) + Nes = 9. (48) 


| 
| 
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Here, Feg¢ is the effective area of the servo system's 
membrane, Nos is the jet reaction of the nozzle R.s, 
defined by the expression 
8 
: where Fes is the effective area of the servo system's 
a 7 nozzle, 
iy f It is easy to obtain the dynamic equation of the 
device from the system of equations (32) and (43) 
4 through (46) if we neglect small terms and bear in 
a remy 7 6 7) 70a atm. mind that, in the square-root extractor, Py = Py: 
Fig. 7. 20,7 
+ [AE Ps — — 
, _%75 ‘41 
ae |S — P,P,. ) 


k k dt jdt 
prs m? Here, the following additional notation was 
used: 


2 The dynamic equation just obtained is essenti- 
ated ally nonlinear. 


ge" In order to improve the dynamic properties of 

the square-root extractor, it is necessary to strive to 
“0b decrease the time constants T, and T;, to decrease 
the impedance ratio 04, and to increase the effective 
relative area k, 


a a a6 48 Axgen® 
Fig, 8. 8. Experimental Results for the 
Multiplier-Divider. 


The multiplier-divider was tested for hysteresis, 
for linearity, for functioning in the multiplying and squaring modes, for carrying out concrete mathematical 
operations and for operating with factors interchanged. 


The results of the hysteresis testing are shown in Fig. 4. The tests were carried out in the following manner. 
With a fixed pressure, Ps = 1 atmosphere, the relationship P, (Py) was read off, with P, acting as a parameter, 
taking the values of 1 and 0.5 atmospheres. The characteristics were read off with P, increasing from 0 to 1 
atmosphere and then decreasing from 1to 0 atmospheres, with a delay at the maximum point for 15 minutes so 
that a residual deformation of the membrane would be produced. The experimentally determined points are 
shown as circles on the characteristic curves, The readings shown were obtained with standard manometers with 
accuracies of 0.2 and a measurement range of 0-1 atmosphere. The experiments established the fact that, for 
practical purposes, hysteresis does not occur. This is explained by the symmetry of the circuits used in the device 
in which, with static equilibrium, the residual deformations of the membrane are mutually equalized, Actually 
used was a flat membrane of membrane material MP. 


Figures 5 and 6 give the results of the tests on linearity in the multiplication mode. These results were 
obtained by reading off the relationship P,(P,) with P, as a parameter, and the relationship P, (P,) with Py as a 
parameter. Two identical families of lines were obtained, with the dispersion of the points not exceeding half 
a scale division on the standard manometer, which comprises an error of 0.35% (together with the error of the 
standard manometer), 


| 
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Figure 7 shows the results of testing the squaring mode, The deviations of the experimental points from 
the parabola are not in excess of half a scale division of the standard manometer, and the device's error in the 
squaring mode is of the same order of magnitude as the error in the multiplication mode. The errors of the 


device in executing the other operations deriving from Eq. (1) were of the same order of magnitude as the experi- 
mental results just discussed. 


The tests showed that the natural time of response of the device was 0.5 seconds, This response time could 
be reduced by decreasing the volume of the inter-valve chambers, Response time is impaired by increasing the 
capacity at the device's output. If such a large capacity is connected then, for attaining a fast response time, 
it is necessary to add a power amplification stage. 


Below are given the results of concrete mathematical operations, on the basis of which one may judge the 
accuracy of the device's performance and its performance with the factors Py and P, interchanged: 


0x0=0.000, 0.3x0.3=—0.060, 
Ox 1=-0.000, 0.2x0.3=-0.060, 


1~0=0.000, (?=0.000, 
1x1=1.000, 12=1.000, 

1x0.5--0.498, 0.72=0.488, 

0.5x1=0,501, 0.5%=0.250. 


It is clear from the examples that the device executes the mathematical operations with three-place 
accuracy, and possesses the property of commutativity of the multiplicative factors. 


9. Experimental Results for the Square-Root Extractor. 


The squace-root extracter was tested for hysteresis, for reproducibility of its indications, for carrying out 
concrete square-root extractions and for response time, The experiments showed that there was virtually no 
hysteresis, that the device's indications were sufficiently stable and were reproduced during lengthy testing. 


The accuracy with which the operations were performed corresponded to an instrument of 0.5 accuracy 
(with measurements made on a standard manometer of 0,2 accuracy), The inherent response time (i.e., the 
response time of the device with no load on the output) was 0.5 seconds. 


Figure 8 gives the family of characteristics read off experimentally when pressure Py varied and pressure 
P, was held as a parameter, 


Below are given examples of concrete square-root extractions executed by the device, from which its 
accurecy may be judged: 

V 0,04=0.2, V 0,81=0.905, 

V0,09=0.295, V1,0=0.996, 

V0,16=0.403, V0,2x0.3=0.246, 

V0,25=0.505, V0,1x0.9=0.303, 

V 0,36=0.606, V0,5x0.7=0.498, 

V0,49=0.705, V0,7x0.9=0.795, 

V0,64=0.806, V0,7x0.7=0.697. 


In the examples, the input and output quantities have the dimensions of atmospheres, It is clear from the 
examples that the square-root extractor performs with three=place accuracy. 


The universal multiplier-divider and square-root extractor which we have developed should find wide 
application in pneumoautomata and in pneumatic analog computing technology. 


These devices may also be employed in hydro-automation, either directly, or in hydraulic variants, pro~ 
vided only that additional measures are taken to avoid the possibility of obliteration of the valve organs. 
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DETERMINATION OF THE EXTERNAL CHARACTERISTIC AND CALCULATION OF 
PARAMETERS OF A HYDRAULIC COMPENSATION REGULATING UNIT 


Vv. M. Dvoretskii 
(Moscow) 


An expression for the external characteristic of a hydraulic compensation re- 
gulating unit is determined, Calculation of the main parameters of the unit is pro- 
posed. Some experimental results are described. 


Pneumatic regulators working on the balance-of- forces principle,for instance, U.M.S, (unified modular 
system instrumentsyhave recently gained wide use in automatic control of production processes, The superior , 
precision and sensitivity of these devices are due to the absence in them of rubbing parts and to the minuteness 
of mechanical movements of regulator components during the operation (of the order of hundredths of a mm), 
The principles used in the design of pneumatic devices could not, until recently, be applied to hydraulic regulators 
which usually incorporate slide valves and jets in their design. Difficulties in controlling small flows of mineral 
oils, precluded the application of above principles, These difficulties were overcome as the results of research 
work [1, 2] conducted in the IAT AN SSSR. This work led to the construction of a hydraulic compensation re- 
gulating unit incorporating "nozzle-vane"-type regulators and controlled-throttle packages. 
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Input APin = 0.4 kg/cm* 


Output__/ X = 0.25 kg/cm? 


L i JO sec 
Fig. 2. 
12 1. Diagram and Principles of Operation. 
‘a va The hydraulic regulating unix enables one to have two types 
| J \e of regulating action: static and isodromic. The diagram of the 
*. | / 7 unit is shown in Fig. 1. 
L 
| Pressure proportional to the instantaneous value of the 
7 a 7 7 regulated quantity, is applied to chamber 04, while pressure Pr, 
2 : 77 7 YA proportional to the required value of the regulated quantity, is 
FA r) applied to chamber 0. The effective areas of membranes My, 
Em VA VA Mg and Mg which form the comparison element are chosen in such 
B. ts | / / a manner that when a pressure difference exists between chambers 
0 and 04, rod Ill, is acted upon by a force whose value and 
I YB - direction are determined by the value and sign of the input signal. 
a The displacement of the rod changes the distance between the 
9 nozzle and the vane, thus regulating pressure P, in chamber 1. 


Input signal (Pyiea5~ Pre kg/ cm* _ Pressure Py controls, by means of the ball amplifier, pressure P, 
in the working chamber 2, Output pressure P,, acting on mem- 
Fig. 3. brane Mg, provides a negative feedback while pressure Ps pro- 
vides a positive one. By changing the hydraulic resistance of the 
throttling range (amplification factor)cock TR, it is possible to vary the relationship between the pressures in the 
positive and negative feedback loops and thereby vary the slope of the static characteristic of the unit. 


An isodromic regulation is achieved by feeding the working fluid through cock AF into the variable-volume 
chamber 4. The value of pressure P, depends on the quantity of liquid which flowed into the variable-volume 
chamber. When the liquid flows through cock AF,pressure Py and, consequently, pressures P; and Ps vary. This 
ensures the variation of the output pressure of the unit at a rate depending on the flow of the working fluid 
through the isodrome cock. With cock AF closed,pressure Py remains constant and a static regulation is thereby 
achieved, 


For higher values of isodrome time it is necessary to ensure small flows of mineral oil through the cock AF. 
The phenomenon of obliteration prevents this being achieved by simply decreasing the port cross sections. In 
the hydraulic regulating unit, therefore, use is made of controlled-throttle packages designed on the principle of 
multiple throttling of the working fluid through apertures of large enough cross sections (see [2]). 


2. Determination of the External Characteristic. 


If the membrane stiffness is neglected, the condition of the balance of forces which act on the rod of the 
comparison element leads to the expression 


APin (F; — F2) = (P2 — Ps) Fj, (1) 


where AP;,, is the pressure drop between the required and measured pressure chambers; Fy and Fy are the 
effective areas of membranes, and F, > F,. 


The flow of the working liquid through the chambers of positive feedback is, according to Bernoulli's 


Q; = 


K, = 100/TR% 
50) 


where Pg is the pressure in the lower isodrome chamber, f; the 
a 1 area of the TR cock orifice, €; the resistance coefficient of 
the TR cock, f, the. port area of the fixed throttle T,, ¢, the 
resistance coefficient of throttle T,. 
Ki FF The pressure drop through throttle T, is determined 
by the expression 
7 
r) i; Substituting (2) in (3) we obtain 
2 
/ AP, = (4) 
fj + 
/ 
iY As Ps s 
Z 
Ps = AP, + P,. (5) 
a 2 J ‘ 6 


Solving jointly Eqs. (5), (4) and (1) we obtain 
Number of turns of the valve ¢ 


Fig. 4. 
APy (Fi — F2) = (1 — P,F, 
Pin (Fi 2) ( oF; (1 + 
whence 
APin (F1 — F 2) | 
Py = (14 + P,. (8) 
aes To determine pressure Py let us consider the variable~ 
{ volume isodrome chamber. The pressure increment in the 
chamber for an infinitesimal time interval is equal 
dP, = ap, (1) 
Ff 
Fig. 5. where Q, is the working fluid flow through the isodrome 


cock (noting that Q, > 0 when the chamber is filling and 
Q: < 0 when it is emptying),Kj the total stiffness of spring Sz and the isodrome chamber bellows, F, the effective 
area of the bellows, The flow through the isodrome cock will be 


where fj is the port area of the isodrome cock, and ¢, the resistance coefficient of the isodrome cock, Hence: 


P,—P,) ... 
dP, =f, A (9) 
Solving Eqs. (6) and (9) jointly and integrating,* we find 


*In integrating it is assumed that APip is not a function of time. 
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$1 At the initial conditions of t = 0 and Py = Pg, which occur 
"be in the regulating unit, the constant of integration C = Py is the 
initial pressure which: exists at the output of the unit at the 
et instant the signal appears. From Eqs. (6) and (10) the general 
expression for the internal characteristic of the hydraulic re- 
as gulating unit is obtained: 


a Fi—F 
ay 


\ 2g (42) 
+ Po. 


K; Fi,—F 


\ 


= 0.05 cm From expression (11) it will be seen that the output 
pressure of the hydraulic regulating unit contains proportional 
and integral components, The integral component, in a general 
Fig. 6. case, is not a linear function of the input signal. However, for 
a laminar flow of the fluid through the isodrome cock, the 
7, min relation between the rate of increase in output pressure and 


9 the value of the input signal becomes linear. [In expression 


(11) the laminar nature of the flow is accounted for by the 

change in ¢, with the variation of the pressure drop in the 

ia isodrome cock, since the form of the expression, characteristic 
YA for the quadrature condition of the fluid flow, is being preserved 

a in its general aspect). 


La 
V4 The experimentally determined relation between the 
Y output pressure and the input signal for TR = 160% and isodrome 
time Tj = 60 sec is shown on the oscillogram of Fig. 2. 


3. Determination of Regulation Characteristics, 


‘ The relation between the adjustments of the unit throttling 
cS @ & & 4 0 Rh, range and isodrome time) and the positions of adjusting controls 
Fig, 7. are taken as the regulation characteristics of the hydraulic 
regulating unit. 


The throttling range adjustment control (the amplification factor of the static part) of the regulating unit 
was made in the shape of a conical throttle valve or a regulating package [2]. 


Henceforth, for the sake of convenience, we shall use the concept of the amplification factor, which can 
be expressed in terms of the throttling range in the form 


2 100 
~ (12) 


where P§ is the proportional component of the output pressure. 
On the basis of expressions (11) and (12) we can write 
F, Celt 
In the case of the conical throttle valve the expression for the amplification factor will take the form: 


(14) 


where Dg is the diameter of the conical throttle valve seat, S the throttle-valve thread pitch, a the angle at the 


apex of the cone, g the number of turns the throttle-valve makes from the fully closed position, d, the diameter 
of the opening of the fixed throttle T,. 


The ratio €,/€, is a function of the pressure at the input of cock TR and of the throttle port area, and is 
determined experimentally by the flow characteristics of throttles TR and [,. Figure 3 shows the relation between 
the output pressure and the input signal with P, = 0 at different amplification factors of the unit working under 
static condition (Pg = Pg = 0), while Fig. 4 shows the regulation characteristics (experimental and calculated for 
€,/, = 1) of the TR cock. Curves 1 correspond to a = 15°, and curves 2 to a = 10°. Dotted lines denote cal- 
culated, while continuous ones experimental curves, 


When a controlled package is adopted the expression for the amplification factor becomes: 


K,= Fi (145-2), (14") 


where dy is the aperture diameter of the package,m the number of series-connected openings, €, the resistance 
coefficient of one opening. 


The adjustment of the isodrome time is achieved by means of a rotary cock with multiple throttling of the 
liquid through openings in diaphragms [2]. 


In order to plot the regulation characteristic of the isodrome cock it is convenient to make use of the time 
required for the statistical part of the output signal to double, (Determination of the doubling time is clear from 
Fig. 5). 


For the case of the isodrome cock Eq. (9) can be written in the form 


P,—P, 


where d; is the diameter of the diaphragm aperture, nj the number of series connected openings, ¢; the resistance 
coefficient of one opening. 


Integrating (15) between t = 0 and t = Tj and between P, = 0 and Py = APinKg (Fig. 5) and solving with 
respect to T; we obtain 
Ki 


The experimentally determined relation between €; and nj at various values of AP = P,~P,, required for 
plotting the regulation characteristic of the isodrome cock, is shown in Fig. 6 while the regulation characteristic 
for d = 0.05 cm and Kj = 6.5 kg/cm in Fig. 7. The calculated and experimental curves are shown by dotted and 
continuous lines respectively. 


SUMMARY 
1. An equation for the extemal characteristic of a hydraulic compensation regulating unit is obtained, 
2. Expressions for regulation characteristics of the unit are determined, With the aid of these expressions 
it is possible to calculate the basic characteristics of throttling elements similar to the ones used in the unit, 
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THE CALCULATION OF TRANSIENT RESPONSES IN 
COORDINATED AUTOMATIC CONTROL SYSTEMS 


N. L. Raikhel’ 
(Briansk) 


A method is presented for constructing the transient responses in coordinated 
linear automatic control systems when control autonomy is disturbed. The transfer 
functions for the controlled quantities are derived, and the concept of partially auto- 
nomous systems is presented. As an example, the curves are constructed for the trans- 
ient response when the load varies on asteamturbine control system using tapped- 
off steam for heat generation. 


Systems for the coordinated automatic control of several quantities occur in the control of steam turbines 
in electricity generation and in commercial uses of steam, as well as in other domains of technology. Most 
advantageous of these systems, from the standpoints of use and design, are the autonomous control systems [1-4] 
in which each of the controlled quantities is controlled independently of the other controlled quantities, i.e., 
without disturbance of their steady states. However,a strict implementation of the conditions for autonomy, 
particularly the conditions of dynamic autonomy [3, 4] is, in practice, quite difficult. These conditions are 
generally implemented with known approximations, It is therefore useful to investigate the processes in coordi- 
nated control systems when the conditions of dynamic autonomy are disturbed (do not hold), The investigation 
of such processes is the aim of this paper. For this, the conditions for static autonomy are assumed to be met. 
Similar types of calculations allow estimates to be made of the influence on the quality of control of the degree 
to which.autonomy is disturbed, With this, the quantitative deviation from the autonomous state is approximated 
by the magnitude of the relative difference of the servomotor times: 

6, — 0; 


where 6; and 6, are servomotor times. 


To solve this problem we use the frequency methods of calculation which are usually employed in treating 
automatic control systems for one variable [5, 6]. 


Since the conditions for dynamic autonomy reduce to the equality of the servomotor times, the method of 
calculation suggested can be used for judging the influence of these times on the quality of the control. 


In this paper, detailed consideration is given to systems for the coordinated control of two quantities. 
Allusion is also made to possible ways of extending and generalizing the results obtained to systems containing 
three or more controlled quantities, and the concept of partially autonomous control systems is introduced. 


1. Equations of Motion and Autonomy Conditions. 


We consider a linear system for the coordinated control of two interdependent variables (quantities) under 
conditions which allow the influence of lags, parameter distribution and friction to be neglected. We also assume 
that there are no crossed connections.* The differential equations for the free oscillations of the system can be 


*On crossed connections, see, for example, [3]. 
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written in the following manner: 


Py (P) = Tia. + Trae, 


= G1, do(p) = Ga (2) 

+ 1) pa = + (3) 
(Oop + 1) = + 


where g and ¢ are the dimensionless controlled quantities, Py(p) and P, (p) 
are operators on the controlled objects, r44,....,T2_ are dimensionless co- 
efficients, either positive or negative, relating the controlled quantities 

and with the controlling acts py and pg, and are the servo- 
motor time constants, in seconds, byy,.....,bgg are dimensionless coefficients, 
either positive or negative, relating the controlling acts py and py, via the 
controller coordinates €, and €, to the controlled quantities, dy(p) and 

d, (p) are controller operators and p = d/ dt. 


Equations (1) are the equations of a controlled object, the processes 
of which depend on two interdependent quantities [1-4], Equations (2) and 
(3) describe the processes in, respectively, the controllers and the servo- 
motors [3, 4]. 


We open the system (for example, from the object to the output) by 


eliminating the intermediate coordinates, 14, 2, $y and €,, from the equations we have obtained: 


= M + Myo%ein, 
= M + 


where Mis = 324s + Qy7 Qi =O 5p + 1, ais = bjs/ ds (p), (i, = 1, 2). 


In an autonomous control system, the coordinates at the output depend only on those input coordinates 
which refer to just one variable and, hence, the equalities 


will be the autonomy conditions. 


M,, = 0, M,, =0 (4) 


Of the possible ways of meeting conditions (4), the meaningful one in practice is the case when the operators 


Q; = Qs, which is possible when 


0; = 6,, (5) 


and the operators d,(p) and d, (p) are numerical constants: 


d, (Pp) = d, (p) = (6) 


where 5, and 53 are the “nonuniformities" of the control of the corresponding quantities. 
Using the additional notation, Nj, = Tiy4ys + Tig4gs, We get, when conditions (5) and (6) hold, that Ny = 0 


and Ng = 0 or, in expanded form: 


+ = 0, + = 0. (7) 


Thus, systems can be autonomous which are controlled by ideal controllers with conditions (5) through (7) 
being met where, as is well known, the first of these is called the dynamic, the latter the static, autonomy con- 


ditions [3, 4). 


In the sequel we shall assume that (6) and (7) are always satisfied and, to estimate the deviation from 


TABLE 
Modulus and Argument of the Function q(z) = z/(z +1) for z = jw 


, sec mod q arg q.deg w, sec” mod q arg q, deg 

0.00 0.000 90.0 1.75 0.868 29.5 

0.05 0.050 87.2 2.00 0.894 26 .6 

0.10 0.099 84.3 2.25 0.914 24.0 

0.15 0.148 81.5 2.50 0.928 21.8 

0.20 0.196 78.7 2.75 0.940 20.0 

0.25 0.243 76.0 3.00 0.949 18.4 

0.30 0.287 73.3 3.50 0.962 16.0 

0.40 0.371 68.2 4.00 0.970 14.0 

0.50 0.447 63.4 4.50 0.976 12.5 

0.60 0.514 59.0 5.00 0.981 11.3 

0.70 0,573 55.0 5.50 0.984 10.3 

0.80 0.625 51.3 6.00 0.986 9.5 

0.90 0.669 48.0 7.00 0.990 8.1 

1.00 0.707 45.0 8.00 0.992 7.4 

1.25 0.781 38.7 9.00 0.994 6.4 

1.50 0.832 33.8 10.0 0.995 5.7 

condition (5), we introduce the quantity 
6. — 0, 
n= (8) 
which we shall call the degree of autonomy distur- 
-10 10 bance. 
@ m case 
We note that, in this 
Fig. 4. Determination of the expression =O; + mé =0,N 
mé 

1+D(w)m; 1is the hodograph of D(w), the + 
lengths of vectors 2 and 3 are mod[1 + D (w)m) where 1 is = rags (i, 8 = 1, 2), 
for m = 1 and 2; 6, and B, are arg (1 + D(w)m) 
for m equal, respectively, 1 and 2; form #1, 2. System Transfer Function with Load 


the linear scale is increased by a factor of m. Variation. 


We now consider the behavior of the control 
system discussed above when the load varies, We assume that the transient response starts from a steady-state 
condition (i.e., the forced response has zero initial condition) and that the forcing function has a step~function 
character, and may be represented in the form 


F(t) =1 (t), 
where 1 (t) is the unit function, 


We shall consider the behavior of a linear system when only one controlled quantity, say, g, is excited. 
Under the conditions discussed in paragraph 1, the equation for the forced response is obtained from the system, 
(1)-(3) where the forcing function, F(t), is appended to the first equation, By eliminating the controller co- 
ordinates, ¢, and €,, we get: 


Py(p) = + + 1 (0, (919 + 1) = + 
G2 = Tortts + + 1) = e191 + (9) 
With this, and by virtue of the assumption that m # 0, there is a transient response, not only for the quantity 


%, but also for the quantity g, which, thanks to (7), will return to zero in the new steady state, while the new 
steady-state value of g will not equal zero if the control system is not astatic, 


Fig. 5. Design for controlling rotational speed 
and pressure in a turbine with steam tapped off 
for heat generation: 1) speed controller, 2) 
pressure controller. 


Fig. 6. Curves of the transient responses 
in the dynamic autonomous system: 1) 
for variation in rotational speed, 2) for 
pressure variation, 


Dy; (2) = (2) + ys (2)], Dar (2) = — Hy (2) 


for g, excited, and 


In what follows we shall employ a second index 
(subscript) for the controlled quantity to denote the 
excited coordinate, We then apply the Laplace transform 
to Eqs. (9), using a superscript bar to denote the trans- 

formed quantity, for example, (y) = 71 
0 

etc. System (9) transforms to the following: 


Pi (y) + trots + 
= tori + 
(Oxy + 1) = + 


We obtain the transfer function by eliminating 
from system (10) the transformed servomotor coordinates, 
Pi, and jig, and then in turn eliminating the transforms 
of the controlled quantities, Gg and $y. In the first 
case we shall have the transfer function 4, corres- 
ponding to $y and, in the second case, the function 
corresponding to Before writing the expressions 
for 4; and jy, we change from the variable y to the 
variable z = Oy, setting 9 = 6, and letting z +1 =Q. 
These substitutions are justified in the construction of 
the curves of the transient response, and are equivalent 
to the replacement of the independent time variable 
t by anew dimensionless variable t* =t/@ [5). 


We introduce the following notation: 


Ti 1 


— Now, (z) = v2, =u, i=1, 2, 


D(z) = (2), 9(2) = 


It is obvious that, with the notation used, Ny = 
= ny and = np. 
After some transformations, we get 


(11) 


(z) = — (z), Doo (z) = (z) + 1H 29 (12) 


for ¢, excited. 
In formulas (11) and (12), we set 


= 


u, (2) v, (2) 


Hs (2) = Ai (2) 


hy Thus, the transient responses in nonautonomous coordinated 


& control systems are made up of two components, one of which 
is the transient response in the autonomous system, the second 
of which arises from the disturbance, or infringement, of auto- 
nomy. 
© 
) —_—_—— The possibility of expressing the transient response by these 
a 0 20 7) wat two components allows an estimate to be made of the influence 


of the degree of disturbance of autonomy on the quality of 
control, for example, estimating overshoot and length of the 
transient response as functions of the quantity m. We note that 


the transfer functions 4, and ,, differ only by a constant 
factor: 


®,,. 


Consequently, the controlled quantities corresponding to 
these transfer functions differ only by the magnitude of their 
deviations from the steady-state values. 


The curves for the transient responses are constructed in 


a well-known way from the transfer functions (cf., for example, 
(5, 7, 8) and the Appendix). 


3. Transfer Function for a System with a Change- 


over Controller. 


By determining the transfer function of a control system 
with a change-over controller, one can show that, in this case, 
the transient response is also made up of two components, one 
of which corresponds to the autonomous system. In the case of 
a change-over regulator controlling the quantity 4, the 


Fig. 7. Curves for the component pro- equations of forced response will be (1) and (3), To the first 
cesses due to disruption of autonomy: of the regulator equations, (2), is added the forcing function 

a, b) for variation in rotational speed, [9], so that, instead of (2), we shall have the following regulator 
c) for pressure variation; 1) m =1, equations: 


2) m = 2, 3)m <4. di pr (t), de(p)e = (18) 


We first eliminate the regulator coordinates, ¢, and 3, from the system of Eqs. (1), (3) and (13) and then, 
employing the Laplace transform and the constraints adopted in paragraph 1, we obtain: 


Py (y) Qu = tata + 


14 
P,(Y) = Tera + a4) 


(Oy + 4) = + + 
(Oxy + 1) He = + + 
As was done in paragraph 2, we derive the expressions for the transfer functions from these equations: 
(2) = (2) — (2)], 5, (2) = (2), 


@;,(2) = £- (2) = 14,2) — (16) 


(15) 


where 


Aj (z) = — uj (2) % (2), 


H(z) = Aj (2) 


In the expressions in (14) through (16), the second sub- 


script corresponds to the regulator which undergoes change- 
over. 


4. Estimating the Stability of Motion. 


By considering p as some complex variable, y, and using 
the substitution @y = z, we can present the characteristic 
equation of system (1)-(3) in the following form: 


By (z) + mB, (z) = 0, (17) 


where the term Bg (z) = Nyy) (QP: ~ corresponds to 
an autonomous coordinated control system for which m = 0, 
and the term By (z) = ~ arises from the 
disturbance of autonomy. 


With Eq. (17) in mind, we could estimate the stability 
of the nonautonomous coordinated control system by means of 


Fig. 8. Curves of the transient responses 


in a nonautonomous system: a) for the Nyquist-Mikhailov criterion [5] which, as applied under 

variation of rotational speed, b) for our conditions, would have to be thus formulated: a non- 

— variation; 1) m = 1, 2) m = 2, autonomous coordinated control system for two controlled 
)m =4, 


quantities will be stable if the corresponding autonomous 
system is stable and if the point (— 1; j0) om the complex z 
plane lies outside of the region bounded by the hodograph of the function 


W (2) m 


as z varies along the imaginary axis (z = jw, ~ © <= w= +), and the characteristic equation of the given system 
is presented in the form of (17). 


The theorem just formulated may be extended to systems with three or more controlled quantities, 


5. On Partially Autonomous Systems. 


With definite relationships subsisting between the servomotor times, the method, presented above, for con- 
structing the transient responses can be extended to coordinated control systems containing an arbitrary number 
of controlled quantities, 


The free oscillations of such a system are described by the following equations: 


n n 
Pi (p) totter (p+ 1) =D) aig 1,2, (18) 
8=1 


s=1 
Without any loss of generality, one can consider that 6;> Op, r = 1,2,.......m— 1. 


Here, the disturbance of the conditions for dynamic autonomy is characterized by n~ 1 numbers, my, 
which, by analogy to (8), we define by the equations 


= (r= 1,2,...,n—14). 


Let there occur a variation of the load, and an excitation of the steady state for the quantity g,. The 
equations for the forced response are identical with Eqs. (18) with the exception of the k'th equation, which 
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has the form: 


Py(P) = Tks + 1 (t). 


8=1 


Proceeding analogously to paragraph (1), we get the transfer functions ;;, defining the transient responses 
in the controlled system. 


We set 
n= (naa 
Hix r=1 r=1 (19) 
Dy (2) + D (2) 


where, for example, C{f) m(*) denotes the product of the polynomial cf) = cH (z) of the complex variable z 
by one of the possible two-factor products, mr), of twoof the numbers m4, Igy...+++eIMg-3- . The number of such 


terms is obviously equal Co ) i.e., the number of combinations of n= 1 things taken 2 at a time. The other 
terms in the numerator of (19) are defined analogously. In the denominator of this expression, 


Do (2) = — (QP 2 — (QPn— Nan), 


n 
where Nj, = ps Tir 4rs- The term D®* (z) is a polynomial with the same structure as the numerator, 


f=] 
We note that, if the conditions of static autonomy hold, Nj, = 0 fori # s. 
It can be shown that the transfer functions havethe following forms: 


Dix (2) = Hix (z) #*) 
Dyx(z) = Hx. (20) 


If the conditions of dynamic autonomy hold, then 


For practical applications, there might be interest inhering in the case when the conditions of dynamic 
autonomy hold for all the controlled quantities but two, for example, g; andg,. In this case, for r 2 2, all the 
numbers, mz = 0, my # 0, and all terms containing at least one of the numbers, Mg, Mgy.+..+,IN,.45 drop out of 
expressions (19) and (20). The transfer functions thus obtained will be structurally identical with the transfer 


functions of a nonautonomous coordinated system for controlling two quantities and, consequently, the method 
presented above for computing the transient responses can be applied to them. 


Partially autonomous systems, in which only some of the numbers, m,, vanish are of practical interest 
since, in many technological processes which depend on several controlled quantities, it suffices to attain the 
independence of only a very few quantities, since the others can have only a second-order effect. 


APPENDIX 


Construction of the Transient Response Curves. 


A frequent case of coordinated control in practice is the case when the controlled object has the same 
load for each of the controlled quantities*, i.e., 


P\(pP)=% p+i1, Pal(p)=—tap +i, 
*For example, a steam turbine with steam tapped off for heat generation. 
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where ry and rz, are time constants. 


In this case, the expressions for wy(z) and wg (z) are of the same type so that, omitting subscripts, one can 
write: 


where r* = 7/0. 


To construct the transient response, we allow z to vary along the imaginary axis, i.e., z = jw [5, 7]. We 
set arg (z +1) (w), arg(r*z+1)=¥%*(w). Then tan ¥ =w andtan¥* = r*w, 


It is easily seen that 


mod = cos ¥, mod = cos 
mod [w (z)] = cos cos 4°, arg[w(z)] = + 9’). 


Thus, the calculation of the modulus and argument of the function w(z) can be carried out with trigono- 
metric tables or with a slide rule. The modulus of the function v(z) differs by a constant factor from the modulus 
ofw(z), and its argument differs by 180°. 


For the determination of the modulus and argument of the functions uy and u, it is convenient to use nomo- 
grams, bearing in mind that mod uy = 1/ mod (1 + yy), arg uy, =— arg (1 + vy), and analogously for u. The nomo- 
gram‘constructed for computing the modulus aid argliment of the expression (1 + v) (subscripts omitted) is the 
union of two polar plots, the distance between the poles of which equals unity (Fig. 1). From the point (1, j0) 
we plot the modulus, v( w), at the angle arg v(w); the segment joining the point (O, O) to the end of this vector 
equals 1 + v = 1/u. On the lines of the polar plot about the point (O, O) are plotted the moduli and arguments 
of the function u(w). If either Ny or Ng is several times greater than unity then, instead of the polar plots of 
Fig. 1, one may use a simple graphical construction (Fig. 2). 


From the point (O, O) we draw segment 1, the vector vy at the angle 8, = arg v. The segment ab (vector 2) 
equals 1 + v = 1/u. Consequently, mod u = 1/| ab| and arg u =~ B3. 


The determination of the modulus and argument of the expression (nyUg + ng) is also carried out graphi- 
cally, as shown in Fig. 3. From point (O, O) we draw segment 1, the vector nyu, at the angle By = arg mus; 
from the end of this vector we draw segment 2, the vector nu, at the angle By = arg ngly. The length of vector 
3, which joins point (O, ©) with the endpoint of vector nguy, equals mod (nyu, + ngly) and angle B = arg (nyu, + 
+ Dgty); here it is also possible to compute D (w) = (myug + Ngly)q, using the table for determining gq. 


In order to determine the modulus and argument of the expression 1 + D(z)m, we construct the hodograph 
of the function D(w), noting on it the points corresponding to the values w = uw, Ww = We, etc, From the point 
O4, one unit to the left of point O (Fig. 4), we draw segments to the points just noted on the hodograph; the 
lengths of these segments, measured in the scale employed, will be the desired moduli, and the angles of these 
segments, measured from the line 0,0, will be the arguments of the function 1 + D(w)m for m =1. For the 
case when m # 1, point O, is displaced to the left of point O, not by unit length, but by the distance 1/m; in 
this case, the scale used in determining the moduli will be increased m- fold. 


By carrying out the preliminary computations just discussed, we find the moduli and arguments of the 
functions Aj (w) and Hj;(w), from which it is possible to find the transfer functions which define separately the 
component of the transient response of the dynamic autonomous system and the component due to the disruption 
of autonomy for some value or other of m. The curves of the transient response can be constructed by the use 
of, for example, the method of trapezoidal characteristics (7, 8). 


Example. Figures 6-8 show the construction, by the method presented here, of the curves for the transient 
responses in a linearized system for controlling a steam turbine with steam tapped off for heat generation (Fig. 5) 
with abrupt variations in load, the following initial data is given: 


= 50sec , = 2.5, T13 = 3.0, 
5sec, = 12.0, = —7.0, 
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4, >= —0,.34 45, 
ay, == —O.58 Ay = +1,21. 


Corresponding to these data we have 
A = —2.06, ny=0.327 nm, 0.073 =0.208 


Figure 6 gives the curves for ci, (t) and (t), the deviations of the controlled parameters, and as 
fractions of the control nonuniformities in the dynamic autonomous system, 


In Figure 7, a,b and c,are constructed the curves Nyy (t), Nyy (t) and hgs (t), corresponding to the transfer 
functions Hy, Hyg and Hgg; these curves represent the deviations of the controlled parameters which are due to 
the disruptions of autonomy, 


The resultant curves for the deviations of the controlled parameters, gq and @gg (Fig. 8, a and b), are con- 
structed in accordance with the formulas; 


pir (t) = (t) + (4), Pan (t) = [Gan (t) + (¢)]. 
The curve for g(t) differs from the curve for hy (t) by the factor ~ eng; the curve for gq; (t) is obtained 
by multiplying by 
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CONTACTLESS RELAY WITH TRANSISTORS 


F. L. Varpakhovskii and R. A. Lipman 


The paper deals with a contactless relay using transistors, designed as a two- 
stage dc amplifier with strong positive feedback due to the total emitter resistance; 
an analysis of the circuit's operation, fundamental design equations, and results of 
experimental verification are given. When P4, P8 and P6 transistors are used, this 
ne 4 Senpiitier produces an output of 100-200 watts with an amplification factor of 
10°-10". 


When a transistor operates under conditions of continuous output current (voltage) variation, the output 
power is limited by the maximum permissible collector dissipation which is governed by thermal considerations. 
In pulse (relay) operation, the transistor’s output circuit has only two fixed critical conditions. In the first con- 
dition the transistor is maximally “closed,” in the second it is maximally “open.” In both cases the power dis- 
sipated in the transistor itself is negligibly small. If the time of transition from one state to another is sufficiently 
small compared to the duration of the stationary condition, then, obviously, the maximum power delivered to 
the load will be limited only by the quantity 


max! c max 


where U., ax {8 the maximum permissible collector voltage and I, »,x is the maximum permissible collector 
current. In this case, the heating-up of the transistor will be negligible and, consequently, the maximum ambient 
temperature can be taken to be the permissible temperature of the p-n junction. In this way, the pulse mode of 
operation permits one to increase considerably (in seyeral cases) the output power of a transistor amplifier. 


9 


Pind 


Fig. 1. Principal schematic of the contactless relay. a) using two transistors, b) 
using three transistors. 
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Fig. 2. Qualitative characteristics of the circuit in Fig. 1. 


The pulse mode of operation of a transistor can be applied not only instrict relay systems but also for 
continuous control (regulation). In the lattercase, a continuous variation of the output power’s mean value 
occurs, due to the continuous variation of the pulses’ mark-space ratio (analogous to thyratron and contact-relay am- 
plification with continuous control), In certain cases (* vibratory” regulators) such a mode of operation is ensured 
by the control system itself, 


In order that an amplifier operate in a pulse mode, it must have a relay-type input-output characteristic 
regardless of the law governing the input signal’s time dependence, The present work examines the circuit of 
a contactless relay, i.e., a transistor amplifier with a relay characteristic, obtained by means of positive feed- 
back. If P4, P3 and P6 transistors are used, such a relay produces an output of 100-200 watts with an amplification 
factor of 10°=10°, 


The circuit of a contactless relay (Fig. 1) is a two-stage dc amplifier with both stages connected to a common 
emitter. The amplifier contains positive current feedback due to the emitter resistance Rg; Rg = Ry is the load 
resistance, and Rj is the internal resistance of the signal source, 


An analysis of the circuit's mode of operation, together with a derivation of the basic design equations, are 
given below. 


Let us examine qualitatively the input-output characteristic of the circuit in Fig. 1,a. Initially, we will 
consider the input quantity to be the voltage u; = e;— ug (designations and the sign convention for currents and 
voltages are shown in Fig. 1). When u, = 0, transistor Vy is “closed” and ugy = 0, while transistor V, is com~ 
pletely “open,” 1.e., 


E—(thep+%) 
0, = = . (1) 
Ucn ies Ry, “i. 
When the transistor is open Uce, + up « E. For conditions (1) to be satisfied, it is necessary that the relation- 
ship 
Ri (2) 
be satisfied, where 
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Fig. 3. Experimental characteristics of PZV transistors, used in the design of 
the circuit in Fig. 1, and the circuit for obtaining these characteristics. 


is the current amplification factor of transistor V, when uc, = 0 and icg = E/R,. When uy; increases, transistor Vy 
opens and current icy increases, Current ij) can be considered to be constant and equal to 


in all of the circuit's operating conditions, since for all conditions ug, + Up « E. 
Hence, with increase in icy current 


iba = — er = — (5) 


As long as condition 


2 E 
tb: > BAL (6) 


is satisfied, and consequently uc, < 0, transistor V, remains open, and the value of the load current i; = icg 
remains practically constant and is determined by Eq. (1). The slope of the input-output characteristic at the 
same time is equal to zero (the region A'B’ of the curve ic, (uj) in Fig. 2,a). 


When, due to increase in uj, inequality (6) breaks down, the current ics begins to decrease while the 
collector voltage of V, changes sign and begins to increase, At the same time, the absolute value of the char- 
acteristic’s slope increases sharply. With further increase in uj,the load current decreases rapidly (region B'C’ in 
Fig. 2,a) until the time when transistor ‘V, is “fully open,” i.e., until the time when the voltage u,y becomes 
zero and then changes sign. 


With further increases in u; the load current icg continues to decrease, but much more slowly, since for 
Ugy < 0 the amplification of transistor Vy drops sharply (region C'D* in Fig. 2,a). It should be pointed out that 
“complete closing” of transistor V, (ue, =,0) cannot be achieved in the given circuit, because the voltage 
Ucei = Yeg does not become zero for any values of uj and Ihy, The larger the value of Ry, the smaller is the 
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Fig. 4. Calculated dependences of the firing and release voltages 

on various parameters of the contactless relay shown in Fig. 1: a) 

E = variable, R; = 110 ohms, Rj = 75 ohms, Ry = 1.7 kchm, 

Rg = 1 ohm; b) Ry = variable, E = 50 v, Rj = 75 ohms, Ry = 1.7 kohm, 
Ro = 1 ohm; c) Rg = variable, E = 50 v, Ry, = 110 ohms, Ry = 75 ohms, 
Ry = 1.7 kohm; d) Rj = variable, E = 50 v, Ry = 110 ohms, Ry = 1.7 
kohm, Rg = 1 ohm. Experimental points are denoted by circles, 


voltage drop ucey on the open transistor V; (Fig. 2,b), Therefore, for effective cut-off of the output transistor the 
value of Ry should be as large as possible, The value of Ry is limited by inequality (2) which must, therefore, 

be satisfied with a minimal margin allowing for possible variations in parameters, The higher the current ampli- 
fication factor By, the more effective is the cut-off of the output transistor, since, in accordance with (2) 


(Ri)max = 


Experiment shows that practically complete cut-off of the output transistor (type Pg and Py) occurs already 
when 6 > 10. In those cases when the gain is insufficient, one can make use of a so-called “composite transistor™ 
[1], as is shown in Fig. 1,b. The application of a “composite transistor” also leads to a considerable increase in 


the relay’s current and power amplification factors, and an increase in the circuit's efficiency due to a decrease 
in the losses in Ry. 


Obviously, complete cut-off of transistor Vz (Ueg = 0) can be achieved by introducing between the collector 
of Vy and the base of Vz (Fig. 1) a voltage ¢; = Ucey, of the order of a few tenths of a volt, The introduction of 


voltage €;, as it will be shown below, permits one also to ensure the circuit's normal operation with high ambient 
temperatures, 


We shall further assume that when transistor Vy is open, practically complete cut-off of the output transistor 
is ensured by one means or another, and, therefore, condition 


Iho == Ico, (7) 


is satisfied, where I, is the relay’s minimum load current and Icg is the collector current with open-circuited 
input of transistor 


In this way, the dependence of icg = ij, on uj = €j — Up is represented by the curve A’B'C"D* (Fig. 2,a). 
The sharp change in the slope of this curve occurs when the voltages cy and Up, change signs. Figure 2,a also 
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Fig. 5. Experimental characteristics of the contactless relay circuit of Fig. 1,a 


(transistors PZV, E = 50 v, Ry, = 110 ohms, Rj = 75 ohms, R, = 1.70 kohm, Ry = 1.0 
ohm) under various ambient temperatures. 


shows the feedback characteristic, i.e., the dependence of the feedback voltage 


Up = Roig = Ro (ier + ie) (8) 


on the load current ij, = icp. The current ieg = icg + ipg monotonously increases with increase in icg, while current 
ies = icg + ipy is maximum when transistor V, is cut off (ip, » 0) and monotonously diminishes with increase in 


As the coefficient 6, increases, the feedback characteristic tends to the dotted straight line OM, deter- 
mined by the equality 


Ug = (9) 


Addition of the abscissas of the curve ic, (u;) and the feedback characteristic yields the characteristic 
determining the relationship of the output signal and the control signal ej = uj + ug. If, in a finite region of this 
characteristic,the following condition is satisfied 


Koco = SRoKy > 1, (10) 


where S =0 ic, /8 uj is the slope of the characteristic relating the input and output signals in the case of no feed- 
back, and RgKy = Rgd ig /Dicg is the slope of the feedback characteristic, then the resultant input-output char- 
acteristic will have a relay shape (curve ABCD, Fig. 2,a). In the regions AB and CD,Kgc < 1; in the region 

BC » Kec > 1, while at the points B and C, where Kg, = 1, the contactless relay “releases” and ° fires." 


Making use of families of input and output characteristics of each transistor, for any given values of cir- 
cuit parameters, it is always possible to obtain graphically the resultant characteristic ABCD. However, such *~ 
graphical plot is comparatively cumbersome and impractical since it requires the knowledge of four sets of 
characteristics,the determination of which is a laborious problem in itself. 


For practical purposes it is desirable to obtain expressions relating the coordinates of the fire and release 


points of the relay IL, I», ef, ej to the circuit parameters E, Ry, Ry, Ro, Rj, making use of as little data con- 
cerning the transistor characteristics as possible. 


The properties of a transistor at low frequencies can be completely described by two families of char- 
acteristics; 


1008 


B= ic /ip=B ic), (11) 
= Ue/ib== (Ue, ic). (12) 


Both functions (11) and (12) (current amplification factor 8 and input resistance r;) depend, to a small 
extent, on Uc (when u, = 0). On the other hand, for the operating conditions of the circuit which are of interest 
to us the quantities ucy and ugg are close to zero: when release occurs Ucar 0, Ucy= Ugg Ug, hasa value of the 
order of tenths of a volt; when the circuit triggers,ucy s 0. At the instant of firing, the voltage uc, is close to 
E; however, at this time transistor V, is cut off and its characteristics in this state are not required for design 


purposes, 


In this way, in determining the coordinates of the firing and release points one can, under the condition 
Uc = 0, ignore the dependence of 6 and rj on collector voltage by assuming that 


 OB/duc=0, (13) 
rj (i, Or, /Ouc=0. (14) 
In this case characteristics (13) and (14) should be determined with ug = 0, This last condition simplifies 
the experimental determination of 8 and r; and the selection of transistors, since for this we can make use of the 


simplest form of circuit (Fig. 3) requiring only one power supply. Moreover, the characteristics can be deter- 
mined for values of current right up to the maximum ratings, with practically no heating of the transistor. 


An approximation of transistor static characteristics, similar to that of (13) and (14),is proposed and widely 
used for design in {1}. In this particular case such an approximation is particularly acceptable because for those 
operating conditions of the transistor which are of interest to us the typical value of collector voltage is fixed 
(Uc 0), 


Making use of (13) and (14) for any operating conditions of the circuit in Fig. 1, we can write 
= (Ri + + (ier + ter) Ro 
i i i 
(Ri + ri) + [ica (4 + Go) + + Ro (15) 


The relay releases at the instant when condition (6) is violated. Then, in accordance with (1) and (5) the 
following relationships hold 


r E r 
Substituting (16) into (15) and neglecting the quantity 1/ 6482 as small compared to unity, we obtain 


At firing, in accordance with (7) and (5), the following expressions are effective 


f E 
Hence, 


f 
f E(Ri+ri) , ER 


Expressions (15a), (15b), (16) and (17) relate the coordinates of the fire and release points to all the circuit 
parameters and transistor characteristics and, therefore, permit one to solve ahy circuit design problem, In 
addition to this it is necessary to have one characteristic Bg(icg) of the output transistor and two characteristics 
By (icy) and rj (icy) of the input transistor. 


*Here and further on, the superscript "r" indicates the value of a variable quantity at release, and the super- 
script "f{ indicates its values at firing. 
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Fig. 6. Experimental characteristics of the contactless relay shown in Fig. 1. 
Rj = 150 ohms, ~~~ R; = 100 ohms, — -— * Rj = 50 ohms. 


Figure 4 shows the calculated and experimental relationships between e] and ef and E, Ry, Ry and Ry of 
the circuit in Fig. 1, a,using PZV transistors. In the calculations using Eqs. (15)-(17), the transistor cRanepacienien 
shown in Fig. 3 were used. 


The condition for the existence of a relay regime can now be expressed in the form 


(18) 


The value of Ae; determines the sensitivity of the relay (any value of ef can be set up by means of a proper 
bias). 


we neglect the dependence of 6, Bz and r; on the collector current, and assume that Bi I bs wei 
Bs Bz By and ri wrj then from (15a), (15b), and (18) we will obtain 


7 (Ro (19) 


Here the condition for the existence of a relay regime will take on the form 


>. (10a) 


It is not difficult to show that under these assumptions (10a) is equivalent to (10), 
From the point of view of increasing the sensitivity of the relay, Koc should be picked as close to unity 


as possible. As usual, this is restricted by the instability of the circuit parameters (mainly the transistor char- 


acteristics). For given variations of the parameters, we can, making use of (15a), (15a) or (10a), design the 
circuit in such a way that the relay regime is preserved under all conditions. 


The power gain of the relay can be defined as 


E*/R 
Ky / (20) 
and can be expressed in the approximate form 


The relationships obtained above are valid under the condition that the Icy of the output transistor is suf- 
ficiently small. This condition is violated at high ambient temperatures, At the same time, the relay's current 
multiplication and the quantity Ae; decrease, 
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However, the main unpleasant effect is the increase in the collector power dissipation when ug = 0. This 
effect can be minimized by an additional negative voltage between the emitter and base of the transistor, thus 
cutting it off [2]. In the case of the input transistor this condition arises in a natural manner. When the control 
signal is equal to zero, transistor V, is cut off by the negative voltage: 


— tp = — -(14-g-). (22) 


To cut off the output transistor, where an increase in collector dissipation is especially undesirable, one 
can introduce a compensating voltage, from a separate source €;> Uce}, between the base of V, and the collector 
of V, (Fig. 1,a). Moreover 


Yes = Ucer — & <0. (23) 


The effect of ambient temperature on the relay*s characteristic in the presence and absence of the com- 
pensating voltage e; is shown in Fig. 5. For an ambient temperature f = 20°C (Fig. 5,a) the introduction of a 
voltage e; = 0.5v hardly changes the relay’s characteristic. The firing point becomes sharper, however, when 
€, = 0;also,the minimum load current is sufficiently small (I, = 1.5 ma, current multiplication is 300), When 
t° = 40°C (Fig. 5,b) and e; = 0,the nominal collector current increases considerably, reaching a value close to 
the firing point of about 15 ma. The introduction of e, = 0,5 v sharply improves the sharpness of the firing point 
and decreases the value of the minimum load current to 1.5 ma. However, the circuit will operate under these 
conditions when e; = 0 if the value of the resistance Rg is increased. 


When the temperature t° = 75°C and e; = 0 the relay regime degenerates and the circuit ceases to operate, 
When t® = 75°C and e; = 0.5 v (Fig. 5,c) triggering is very sharp and the minimum current is 2,2 ma. 


In this way, the introduction of a constant voltage €; = 0.5 v guarantees the circuit's completely normal 
operation at comparatively high ambient temperatures, 


Variations in the temperature lead only to a displacement of the “loop” in the relay’s characteristic parallel 
to the abscissa (for negligible variation in Ae;j). In many cases this displacement is insignificant and, if necessary, 
can be compensated by introducing a temperature-dependent bias. What is significant is the fact that the maxi- 
mum output power of the amplifier does not decrease with increase in ambient temperature. 


In the case of operation into a resistive-inductive load, the collector voltage can, at the instant of release, 
be considerably higher than the supply voltage. If, at the same time, the value of the supply voltage is chosen 
close to Uc maxe the transistor can be damaged. To eliminate overload,the resistive- inductive load should be 
shunted by a diode (or a capacitor) as shown dotted in Fig. 1,b). In this case the inductance “discharges” through 
the diode; there is no overload, and, consequently, the maximum permissible peak power of the transistor can 
be fully utilized (E = Uc max: E/R, = Ig max)> 

Figure 6 shows the experimental characteristics of a relay built according to Figs. 1,a and b using P4, P3, 
P6 and P8 transistors. The ‘output power of the relay shown with its characteristics in Fig. 6,b is 200 watts 
(5a, 40 v). The power dissipated in the output transistor in the “open” state is equal to approximately 3 watts 
(Ucez = 0.6 v, Icg = 5a) while in the “closed” state itis 0.14 watt (Uge, = 40 v, Icg = 3.5 ma). Current li- 
fication is 1400, The relay’s sensitivity for a source resistance Ry = 150 ohms is Ae; = 0.33 v or AP; = Aej / ARiw 
ss 1.85x 10~ watts. The power gain in this case is K, # 1.1 x 10°, 
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DETERMINATION OF HEATING OF ELECTROMAGNETIC RELAY WINDINGS 


M. I. Vitenberg 


The heating of electromagnetic relay windings is investigated. It is shown that 
the generally accepted heat-emission coefficients of relay windings and equipment are 
unfit for the design of small relays. A curve is plotted and formulas are given for 
determining the dependence of the mean coefficient of heat emission on the size of 
the winding’s cooling surface. 


The concept of the specific rise of the winding temperature, characterizing 
the relay’s construction from the point of view of its heat emission, is introduced, 
and approximate formulas are given for determining the rise in winding temperature 
and the heating time constant as functions of cooling surface, weight, and dimensions 
of the relay. 


Under steady state thermal conditions the heat produced by current flowing in a relay winding is dissipated 


in the surrounding medium. Emission of heat by the coil occurs by means of radiation, conduction, and con- 
vection processes. Overheating of relay windings is usually quite small and, consequently, for a practical design 
one can make use of a mean coefficient of heat emission (mean heat transfer coefficient), The outer and 
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Fig. 1. Curves showing the mean rise in the temperature of 
relay winding as a function of supplied power for 6, = 20°C: 
1) type RSM without cover, 1") the same with cover, 2) type 
RMU, 3) type RS-13 without cover, 3") same with cover, 4) 

type RPN, 5) type RKN, 6) type KDR. 
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Fig. 2, Curves relating the mean coefficient of heat emission to the 
coil diameter when 9 = 50°C and 64 = + 20°C, 


inner surfaces of the winding are not cooled in the same way. The outer cylindrical surface of the winding is in 
the path of a rising air flow and is well cooled, while the inner (cylindrical or rectangular) surface of the winding 
gives off its heat to the steel core of the relay through a layer of insulation (the bobbin) separating the winding 
from the core, The cooling of the inner surface of the coil depends on the thickness and thermal conductivity 

of the layer of insulation (or spool), the size of the air gaps between the winding and the core, and the thermal 
conductivity of the relay’s magnetic circuit. The coil’s end faces are usually small compared to its outside 
(cylindrical) surface; moreover, they are, in general, covered thickly with paraffined paper, (or plastic) so cooling 
through the end faces can be neglected. However, in those cases where the coil diameter is greater than its 
length, the cooling end surfaces must be taken into account. 


Temperature variations in a relay winding made of enamel-insulated wire are small. The temperature of 
the outer surface of the winding in a RKN-type relay is usually 2-5% below the winding’s mean temperature, 
The maximum temperature of the relay winding is approximately 4-8% above its mean temperature, 


The mean rise of the temperature in the winding of a directcurrent relay above the temperature of the 


surrounding medium, under steady-state thermal conditions, can be expressed by means of the following approxi- 
mate equation [1, 4]: 


P , 
(a) 


where P is the power consumed in the coil, in watts, Sp is the area of the outside cooling surface of the coil, in 
cm’, S; is the area of the inside cooling surface of the coil, in cm”, a is a coefficient depending on the heat 


given off through the relay core, and q is the coefficient of heat emission (heat transfer coefficient) of the coil, 
in watts/cm*T. 


Experiments have shown that for a sheathless (wound directly on the core) coil the value of the coefficient 
a in an RKN type relay is usually somewhere within the limits of 0,9-1.3, In a large relay with a sheathless 
compounded winding, the coefficient may be as high as 2,4[1}. The size of the coil’s internal cooling surface is 
usually considerably smaller than the external one, for which reason variations in the value of the coefficient 
a have a comparatively small effect on the final result of a design. 


In practical engineering calculations of the heating of direct-current electromagnetic relays, the design 
size of the coil’s cooling surface is usually taken to be equal to the sum of the coil’s external and internal cooling 
surfaces, In this case,the rise in the temperature of the relay's winding is found by Newton's formula 


(2) 


where qy is the mean (conventional) coefficient of heat emission of the relay’s coil, and Sg is the coil's design 
cooling surface, 
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Fig. 3. Curves showing the mean coefficient of heat emission of relay 


windings as a function of the coil’s calculated cooling surface, for 
= 50°C and = + 20°C. 


Mean Coefficient of Heat Emission of Relay Coils 


The value of the mean (conventional) coefficient of heat emission of electromagnetic relay coils and 
direct-current equipment is not constant and depends on many factors. The coil's coefficient of heat emission 
depends on the type, form, and color of the coils external surface, the diameter and quality of the wire, the 
material and thickness of the coil"s external and internal insulating layers, the thermal conductivity of the relay's 
magnetic circuit, the temperature of the winding, the relay’s position in space, the speed of circulation of the 
surrounding air, etc, Therefore, an exact value of the mean heat emission coefficient for every particular case 
can be determined only by experimental means. 


The value of the mean coefficient of heat emission for rer’ bay coils and equipment which is cited by most 
authors is somewhere between 0.9x 1078 and 1.4 x 10°® watt/cm 


Thus, for example, E. Iasse [2] points out that the value of the ooeticion of heat emission for coils is 
practically constant; for a dull surface it is equal to 1.2-1.4- -107* watt/cem**c, and only under particularly 
unfavorable conditions falls below this limit. M. A. Babikov [3] acnuies that the coefficient of thermal 
emission for windings with paper insulation is equal to 1,0-1.25- 107° watts/cm**C, F. A. Stupel’ [4] and G, 
Roters [5], give various values for the mean coefficient of heat emission, ranging between 0.94 x 107 3 and 
1.34 x 10°* watts/cm*°C depending on the temperature of the winding (20°C to 80°C) and the quality of the 
thermal contact between the winding and the relay core. In reference [6] the cosffictent mee heat emission of 
the coil'’s outer surface is given values ranging from 0,77 x 107* to 1.34 x 10°* watts/cm**c, 


For practical calculations of the heating in the coils of direct-current relays having a cooling surface of 
approximately 50 to 500 cm’, the values of the coefficient of heat emission cited above apparently gave more or 
less acceptable results. Consequently, in the calculation of the thermal properties of relay windings and equip~ 
ment,the value of the mean coefficient of heat emission has up to now been assumed to be between 0.9 x 107° 
and 1,4 x watts/cm**c. 


However, experience shows that these values of the mean coefficients of heat emission are not applicable 
in the calculation of the thermal properties of small relays because they oe to very arg errors. For example, 
if one takes the value of the mean coefficient of heat emission as 1.2 x 107* watts/cm?*C, then the calceleted 
rise in the temperature of the winding in an RSM type miniature relay, having a cooling surface of 6.5 cm’, 
exceeds the experimentally determined value by a factor of approximately 2.3. Therefore, the value of the 
coefficient of heat emission of small relays is appreciably larger than for normal size relays. 


To investigate this problem, the mean rise of the winding temperature was measured as a function of the 
supplied power on a large number of different electromagnetic relays. Tests were carried out on about one 
hundred samples of electromagnetic relays with twenty different types and sizes of external rotary armatures, 
five prototypes of various sizes of electromagnetic relays, and three separate, large, coreless coils (solenoids). 
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Fig. 4. Curves relating the mean coefficient of heat emission of 
a relay winding to the size of the coil’s external cooling surface, 


The sizes of the cooling surfaces of the coils in the relays tested varied between 1,1 and 240 cm", while for the 
separate coils they were between 1500 and 5000 cm’, 


The relays were freely suspended in still air by wires to eliminate the effect of withdrawal of heat through 
the oven. The mean rise in the temperature of the relay winding was determined by measuring changes in the 
coil's resistance, which was done by means of a Wheatstone bridge. Measurements were made with an ambient 
temperature 09 of about 20°C, 


Figure 1 shows the curves relating the mean pise in winding temperature (above ambient) to the final power 
in the winding for six different types of relays. Over the temperature range 0-70°C these curves have a relatively 
slight curvature. If one takes into account the fact that the rise in winding temperature for different samples of 
the same type of relay usually fluctuates over a range of + 10%, then the curves shown in Fig. 1 can, with 
sufficient accuracy for practical calculations, be replaced by straight lines passing through the origin and the 
points on the curves corresponding to 50°C. These straight lines are shown dotted in Fig. 1. 


Hence, the value of the mean coefficient of heat emission of relay windings can be considered to be 
practically independent of temperature over the range of 20 to 90°C. 


Analogous curves of the winding-temperature increase were plotted for all the tested relays, These curves 
and Eq, (2) were used to determine the mean coefficient of heat emission for each type of relay separately, The 
ratio of the relay's length to its outside diameter for all the relays was between 1 and 3. 


Figure 2 shows the mean curve of the coefficients of heat emission q, of the winding in an electromagnetic 
relay with an external rotary armature, as a function of the coil's external diameter D, plotted from experimentally 
determined results. From this curve it follows that the value of the mean coefficient of heat emission of relay 
coils varies over a wide range depending on the coil diameter. 


Over the range 0.3 < D < 5.0 cm a logarithmic plot of this curve is practically a straight line, Hence, the 
relationship between the mean coefficient of heat emission and the coil diameter can be expressed, over the 
range indicated above, by the following equation: 


q, =aD-, (3) 


where a and b are constants. 


Coefficient a is equal to the mean coefficient of heat emission of a coil of 1 cm diameter. Coefficient 
b is equal to the tangent of the angle of inclination of the linear region of the curve q, = f (D) to the abscissa, 


Calculating the coefficients a and b by the method of least squares, we find ay = 2,98 x 10°? watts/cm**c 
b =— 0,664. Equation (3) takes on the form: 


4 


= 2.98 x (4) 


1015 


| | 
% whom? % 


% when? From Fig. 2 we see that deviations of the value 
> of qy from the mean curve is, for many types of 
PSSh. a relays, appreciable and reaches + 30% Such large 
ai Se y deviations are explained by the fact that heat is with- 
_ SS fe drawn not only through the external but also through 
2 +S 4 ~ coil surface, while the ratios 1 /D and D/d are not 
constant (2 is the coil length, d is the coil's internal 
“> TS Figure 3 shows curves relating qy to the design 
: cooling surface of the coil S,, The dotted curves 


$ 29 570 20H 50110 2008,,cr* show the limits between which the value of qy for 


Fig. 5. Curves of the mean coefficient of heat different types of relays fluctuates. The crosses are 
emission as a function of the size of the total the values of q, for relays protected by covers. 
cooling surface of the relay. Deviations in the value of qy from the mean value 


for most types of relays tested did not exceed + 15% 
The value of the coefficient of heat emission of relay coils should, therefore, be expressed as a function of the 
design cooling surface and not of the coil's outside diameter. Of course, in those cases when the ratio of the coil 
length to its diameter falls outside the range of 1-3, and the construction of the relay differs appreciably from 
the usual construction of a relay with an external rotary armature, the value of q, can differ from the mean value 
by more than +15% 


The curve shown in Fig. Ly has two linear regions with different angles of inctinetion to the abscissa: one 
between the limits of 1-80 cm* and the other between the limits of 200-5000 cm*, The cosfliciens 4 and b 
for the first and second region on the characteristic are mquery se = 6.5 x 10-* watts/cm**C, by = 0.315 
and a, = 3.1 x 107° watts/cem**C, b, = 0.19, 


Thus, the dependence of the mean coefficient of heat emission of a relay winding on the cooling surface 
in the range 1 < S, < 100 cm® can be expressed by the equation 


= 5.5 X 1078 Sy Cy, (5) 
V Se 


while in the range 100 < S, < 5000 cm’, hy the equation 
q's = 3.1103 w (5a) 


In view of the fact that the concept of the design size of the cooling surface is conditional, let us examine 
two extreme cases; 


a) The design cooling surface is equal to the outside cylindrical surface of the relay’s coil Sp. In this case 
the coil’s mean coefficient of heat emission, which we will denote by qg, will be approximately 32 % greater than 
q if Sp = 30 cm*, The curve showing q, as a function of Sp is given in Fig. 4. 


b) The total cooling surface of the relay S, is equal to the sum of the external surfaces of the coil and the 
relay’s magnetic circuit, all in free contact with the surrounding air. Also, the thermal resistance of the coil*s 
internal insulating layer and of the magnetic circuit is assumed to be zero. In this case the relay’s mean coef- 
ficient of heat emission, which we shall denote by qg, will be, for S, = 30 cm*, approximately 17% smaller than 
ay (Fig. 5). 


Comparing and analyzing the curves in Figs. 3, 4, and 5 we come to the conclusion that the most correct 
interpretation of the design cooling surface is the sum of the external and internal cylindrical surfaces of the 
relay coil. 


The equations and curves given above for determining the coefficient of heat emission of electromagnetic 
relay windings can also be used for calculating the temperature rise in the windings of many other analogous 
electromagnetic devices and instruments. 
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Temperature Rise in Relay Windings 


As shown above, in the range 0-70°C the curves of the mean steady-state increase in the temperature of a 
relay’s winding as a function of the power can be replaced, with sufficient accuracy for practical calculations, 
by straight lines passing through the origin and the point corresponding to 50°C. Hence, for each type of relay, 
the ratio of the mean rise in winding temperature to the consumed power, i.e., the mean steady-state rise in 
winding temperature per watt (specific rise in winding temperature) over the given temperature range can be 
assumed‘to be a constant, for all practical purposes. 


The specific rise in the temperature of the relay's winding, in °C/ watt, is 


Knowing the specific rise in the winding temperature of a given type of relay, one can find the mean 
steady-state rise in the winding temperature of this relay, for a given power, in the range 0-70°C: 


= (6a) 


If the winding fills all the winding space on a bobbin, the quantity 9 is practically constant for every type 
of relay; it varies over a small range depending on the diameter and type of wire, the material and thickness of 
the insulating layer between the winding and core, and the thermal conductivity of the coil's external coating 
material, Hence, for a completely filled bobbin, 9, is the thermal characteristic of a given type of relay, char 
acterizing the heat emission of the latter when under steady-state conditions, 


When all the winding space is filled in different types of relays, the quantity 94 is not the same and depends 
on the size of the coil’s cooling surface. Figure 6 shows a curve relating the specific rise in winding temperature 
to the quantity S,, plotted from experimental data obtained when @, = 20°C. The dotted lines are the limits 
between which the quantity 9, can vary from one type of relay to another having a different construction. Over 
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Fig. 9. Curve relating the mean coefficient of heat emission to the 
relay’s weight for W/V =1.5— 2, 09 = + 20°C, 


the range 1 < S_ < 200 cm’ the relationship between 9, and the value of the design cooling surface is practically 
linear; 


9, = 182 
V Se 


The mean rise of the relay’s winding temperature, under steady state conditions, over the indicated range 
of Sc, is 


(7) 


170P 


The slight curvature of 9 = f (P) (Fig. 1) shows that the coefficient of heat emission and the specific rise 
in the relay’s winding temperature depend on the temperature, 


(8) 


From these curves one can find the indicated dependence for any type of relay. An analysis of the curves 
shows that in the range 30 < $< 90°C the specific rise in winding temperature decreases while the coefficient 
of heat emission increases by approximately 2-3% for every 10°C rise in the relay’s winding temperature, 


Dependence of the Rise in Winding Temperature on the Weight and Volume of the 


Relay 


In the design of new types of electromagnetic relays, as well as for determining heating in existing types, 
it is very important to know (approximately at least) the dependence of the rise in winding temperature and the 
heating time constant of a relay on its weight or dimensions, Of course, such a relationship will vary from one 
relay to another, depending on its construction; therefore,measurements were made on relays more or less similar 
in construction. 


Figures 7 and 8 show the experimental curves of the specific rise in the winding temperature of an electro- 
magnetic relay with an external rotary armature as a function of its weight and volume for W/V = 1,5— 2 and 
for 6, = 20°C, The ratio of the core length to its diameter in all the relays tested was between 3 and 8, while 
the ratio of the relay*s weight to its volume was between 1.5 and 2,0, From Figs. 7 and 8 it follows that logarithmic 
plots of the specific rise in winding temperature as a function of the weight and volume of a relay are almost 
straight lines, The dependence of the rise in her ys temperature on the relay’s weight in the range 1< W< 2000g, 
and on its volume in the range 1< V < 1000 cm’, can be expressed by means of the following approximate 
formulas: 


(9) 


g,, watts/cm?°C 
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Fig. 10. Curves of the temperature rise in the winding of an RKN type relay as a 
function of time: 1) U = const, 2) I = const, 
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Fig. 11. Curves of the heating time constant of a winding as a 
function of the relay’s weight: 1)U =const, 2) I = const. 
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where W is the relay's weight in grams and V is its volume in cm’, 


A curve of the mean coefficient of the winding’s heat emission as a function of the relay's weight is also 
shown in Fig. 9, 


In the range 1 < W< 500 g the quantity qy is expressed by means of the following approximate formula: 


Vw 


Heating Time Constant of the Relay Winding 


The heating time constant of a relay winding depends on the heating time, since the relay represents an 
inhomogeneous object consisting of different materials (copper, enamel, paper, varnished fabric, steel, etc) with 
different specific heat capacities and different thermal conductivities. Moreover, if the winding temperature 
rises at constant voltage,the power decreases and the temperature rise slows down, Consequently, an accurate 
determination of the heating time constant of a relay winding represents a very complex problem [4, 5). 


An approximate value for the heating time constant of a winding can be determined from an experimental 
curve of the growth of the temperature of a relay winding. Figure 10 shows experimentally obtained curves 
expressing the time dependence of the increase in the winding temperature of an RKN-type relay for constant 
applied voltage U and constant current I in the winding, as well as a curve of the cooling of a relay winding. 


Exponential curves are shown dotted in Fig. 10, 


From Fig. 10 it follows that the deviation of the curves of the heating of a relay wincing from an exponential 
is comparatively small. The curve of the increase in the winding temperature of a relay at first rises more rapidly, 
and after a time t = r more slowly than an exponential. For a constant applied voltage (curye 1) the heating 
time constant is about 1.5 times smaller than it is in the case of a constant current in the winding (curve 2), 


The cooling time constant of a relay winding is practically equal to its heating time constant when I = const, 
Analogous means were used to obtain approximate values of the heating time constant for relays of different size 
and weight, having cores with length-to-diameter ratios between 3 and 8, The data was used to plot in Fig. 11 
the curves of the dependence of the heating time constant of the winding of a relay with an external rotary 
armature on the relay's weight when U = const (curve 1) and I = const (curve 2). The logarithmic plots of these 
relationships are straight lines and can be determined from the following approximate formulas: 


t, =~ 0.78 VG when U = const, (12) 
when J = const, 


where r is the relay*s heating time constant in minutes and Wis the relay’s weight in grams. 


Solving simultaneously Eqs, (8) and (9) in terms of Wand making use of (12), we obtain approximate 
formulas for the heating time constant as a function of the cooling surface of the relay’s winding: 


3 3 
tu $2 and 1.24/82. (13) 


Examples 


Let us determine the mean increase in the winding temperature of a relay having a coil with a cooling 
surface of 14 cm’, five minutes after switch-on and under continuous operation in an ambient temperature of 
40°C. The resistance of the winding at 20°C is 255 ohms and the supply voltage is 24 v. 


Let us determine the power consumed at the instant the relay is switched on, at a temperature of 20°C; 
P =U*/r = 2.26 watts. 


The temperature increase of the relay winding for a constant power of 2.26 watts is 


1021 


170P 
$= ; =66,3°C, 


V 
The steady-state increase in the winding temperature for a constant supply voltage and an ambient 
temperature of + 20°C [7] (p, 200-202) is 


( 


49 
i+ 334.5 = 54,2°C. 


The increase of the winding temperature of the relay when P = 2,26 watts and the ambient temperature is 
+40°C is 


= 


The steady-state increase in the winding temperature of the relay for a voltage of 24 v and an ambient 
temperature of +40°C is 


234.5 + 0, 49 


The approximate value of the heating-time constant of the relay winding is 


= 0.81 $2 = 4.7 min. 

The increase in the relay's winding temperature five minutes after switch-on, if one assumes that the 

initial portion of the heat-up curve coincides with an exponential curve, is 
= * j= 40.2°C. 
. The power necessary for the relay's reliable operation is 2 watts. Let us determine the approximate weight 

and volume of the relay if the maximum permissible increase in the winding temperature is 60°C. 

From Eq, (9) we find the smallest possible weight of the relay: 


176P\2 
The smallest possible volume of the relay we determine from Eq. (10): 


134P 2 
V = 20 ea’. 
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THE PROBABILITY OF DEFECTIVE ELEMENTS IN 
AUTOMATIC CONTROL SYSTEMS 


S. M. Kuznetsov 


(Moscow) 


Probabilities are determined for defects in, and useful life of, elements in auto- 
matic control systems which are subject to complex stationary disturbing factors. 


Due to the instability of technological processes and to the greatlydiverse conditions of use, the existing 
methods of experimental estimation of element reliability (under standard conditions of use) do not provide the 
necessary accuracy of estimation for elements which are to function in different blocks and links of various 
apparatus, and do not provide a reliable prognosis of the reliability of apparatus which are manufactured in 
successive batches and which are also subject to storage under differing warehouse conditions.* These methods 
entail large expenditures of time and a deterioration of elements and apparatus during the testing process. 


Estimates of reliability based on accelerated testing methods (by forcing the action of the disturbing factors), 
due to insufficient study and lack of theoretical foundations for the test regimen, also do not attain the necessary 
accuracy in estimating the reliability of apparatus under various conditions of storage and use, and entail deter- 
ioration and breakdown of the tested apparatus, which is no less significant than in the former methods. 


In these circumstances, the development of analytica) methods for calculating the reliability of apparatus 
and of individual elements is an important and pressing need, 


Probability of Element Defects.** 


If the functioning of an element is characterized by several independent parameters By, By,......++.By, given 
quantities which characterize the electrical, mechanical and thermal reliability of the materials, the admissible 
values of which are Bj adm, then the proper functioning of such an element may be characterized by the time 
Ty during which the following conditions hold: 


AB AB,(T), 
AB, adm= AB, (T), 


(7), 


where the AB; (T) are continuous functions which characterize the probable variations of the element's parameters 
under the influence of various disturbing factors arising under the actual conditions of manufacture and use, and 
AB; adm are the admissible deviations of these parameters, 


We now determine the relationship between the probability of improper functioning (a defect) of an ele- 
ment as regards one of its parameters and the total probability of element malfunction as regards all of its con- 
stituent parameters. 


* With stringent requirements on element reliability, estimation errors under these conditions might exceed the 
measured quantity by an order of magnitude or more, 

* *By an “element” is to be understood a nondecomposable union of interacting components and materials designed 
to perform the simplest transformations of physical quantities (condensers, resistors, vacuum tubes, etc). 


1023 


We denote the probability of element malfunction as regards each of its constituent parameters by Py, 


Each of these probabilities is a function of T of the form: 


9, (7) aT, 


where the 9;(T) are the differential distribution functions of element malfunction as regards parameter Bj. 


The event that any parameter remains within its admissible limits is the complement of the event that the 


parameter goes beyond these limits. Therefore, the probability of proper element functioning as regards each of 
the constituent parameters during time T will equal Ri = 1— Pj. 


The probability that the element will function properly with respect to all its parameters is, according to 
the law of multiplication of probabilities, 


(4 — Pa) = RyRy... 


The probability of improper element functioning with respect to at least one of its parameters is deter- 
mined from the equation 


Py =1—(1— — P,)...(14— Pa) =1— AR... Ra. (1) 


Since a deviation beyond the admissible limits by at least one of the constituent parameters would disrupt 


the operation of the element, the quantity P, represents the total probability of a malfunction (defect) of the 
element. 


Formula (1) relates the total probability of element malfunction with respect to all the constituent para- 
meters to the probability of malfunction with respect to each of these parameters. For small values of the Pj, 
formula (1) can be simplified and given in a form more convenient for practical computations, By multiplying 
out and then neglecting second-order terms, we will have 


n 
Po= Py +Po+...4+Pr= >) Pi. (2) 
im1 
In the particular case when all the probabilities of malfunction with respect to each of the constituent 
parameters are identical and nP < 1, formula (2) takes the form Py = nP and shows that, for identical probabilities 
of malfunction with respect to each of the element's constituent parameters, the total probability of a defective 
element increases proportionally to the number of its parameters. 


Certain Regularities in the Variations of Element Parameters Under the Influence 


of Disturbing Factors. 


As shown by many investigations [3-10], various disturbing factors arising in the processes of manufacture 
and use of elements can influence not only the initial deviations of the element parameters, but also the speed 
of their variation, giving rise both to reversible and irreversible processes of parameter variation. 


The variation of an element parameter may, in the general case, be represented by a function of the form 


T 
B,(T) = By + {vs aT, (3) 


where Big is the initial value of the element parameter at the moment of time Ty and v; = dB;/dT is the speed 
of variation of the element parameter. 


If the disturbing factors are stationary, the parameter variation of the properly functioning element during 
time Tp is a function of time with continuous n'th-order derivatives. 


1024 


Mexdt) In this case, the function By(T) of any parameter 
on the bounded interval of time, T < Tp, can be approxi- 
mated (with an error not in excess of the remaining terms 
of the expansion) by an n'th-degree polynomial 


= ag + a,T + a,T* +... 4 (8a) 
where the values of the coefficients a9, 44,.......p, a8 is 


known, may be expressed in terms of the derivatives of 
this function 


Poe pees ee 


’ i 
=/(0), =/'(0), a2 = = — J"), 
Fig. 1. expressing, respectively, the initial values of the para- 
meters, their speed of variation, their velocities, etc at 
A(T) Ss the initial moment of time Ty. 


On a bounded segment of time T we may, with 
sufficient practical accuracy, limit ourselves to the first 
two (or three) terms of expansion (3a): 


= Biot +, (Bb) 


4 
a where Vjg and €j9are, respectively, the initial 
deviation, speed and acceleration of variation of the ele- 
Fig. 2. ment parameter, 


Determination of the Components of the Resulting Deviations of the Element 


Parameters, 


The character of element parameter variation in the general case is determined by the physical nature 
(chemical constitution and structure) of the materials, and depends on the physical design, on the technology of 
element manufacture, on a large number of usage factors (parameter regimens, electrical, mechanical and 
thermal loads), and on parameters of the external medium (humidity, pressure, temperature, etc), 


If the disturbing factors essentially affect element parameter variation then, in determining the resulting 
deviations, it is necessary to take into account the effect of the disturbances, not only on the variation of the 
initial deviation, but also on the variations of speed and acceleration of the parameters of the element. 


When there are a large number of disturbing factors, it is necessary to consider the quantities B, v and € 
as functions of many variables; 


B=,(A,, Ag, Am), 
v = f2(A;, Ag, Am), (4) 

= f3(A,, Ag, Am), 
where Ax, in the general case, are disturbing factors which characterize the actual conditions of manufacture 
and usage of the element. Each of these factors may be a constant or a random variable, or a regular or a random 
function of time. 


The variations of these disturbing factors, given the random character of their manifestations, can be given 
statistically by the probability characteristics of the deviations of the disturbances from their nominal (computed) 
values, 


Estimating element reliability when the disturbing factors are nonstationary is very complicated and 
requires special considerations. For the purpose of simplifying the solution of the problem we shall, in the sequel, 
consider variations of element parameters when the disturbing factors are stationary. 
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We now determine the increments of each of the 
functions expressed by Eqs. (4) when the disturbing factors 
deviate from their nominal values by the quantities AA,,. 
If the functions B(A,), v(A,) and € (A;,) have continuous 
n'th-order derivatives, we expand each of them in a 
Taylor series, limit ourselves to the first-order terms (in 

r view of the relatively small argument variations) and, 
after some transformations, we get:* 


= 


AB = a, AA, + + amAAm, 

Av = B,AA, + ©) 

As = y,AA, + ... + 
where AB, Av and Ae are the changes in the components 
of the resulting element parameter deviations under the 
action of the given complex disturbance; a, = 0 f,/dA,, 
By = 0%/0A,, and y, = Ofs/ dA, are coefficients which 
characterize the degree of influence of each of the 

r disturbing factors on the initial values, speeds and ac- 

celerations of the element parameters. 


Within a limited interval, the coefficients of the 
disturbing factors a}, 6), and y, may be considered as 
constants, 


The magnitude and sign of each of these coefficients depend on the character of the element transformation 
by the disturbing act, and may be determined experimentally or by means of numerical solutions of the equations 
which describe the elements" behavior. 


Resulting Element Parameter Deviations Under the Influence of Random Disturbances 


We now deterraine the resulting deviation of the element parameters in the general case when the disturbing 
factors simultaneously affect both the initial deviation and its speed and acceleration. 


If the values of the element parameters at time T are expressed by relationships of the form given in (3b) 
then, if the terms in the right member vary by the quantities AB, Av and Ae, the resulting parameter deviations 
relative to their nominal (computed) values, at any moment of time will equal: 


AB(T) = AB, + (vo + Av) T +(e) + As) T2/2. (6) 
If we substitute the values of AB, Av and Ae in (6) from (5), we will have 
AB(T) = + (v0 + + (co + 72/2, 


or, what amounts to the same thing, 
AB(T) = + (An + Ade) + 1) 


*When element reliability is estimated by forcing the action of the disturbing factors and using a Taylor series 
expansion of each component of the resulting element parameter deviations, it is necessary to take the following 
terms of the expansion into account: 


AB, = 2a, AA, + A? 
Av = 53, 4A, + 53,44}, 
Ae = Ly, AA, + Zy, 44}, 
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where Ay is the nominal (computed) valve of the disturbing factors. The first term on the right side of expression 
(7) characterizes the reversible processes of parameter variation under the action of the given complex disturbance. 
The factor T in the other terms of the expression denotes the presence of irreversible processes of parameter 
variation which are independent of the physical nature of the disturbing factors. 


The heterogeneity and nonconstancy of the disturbing factors and their differing effects on parameter 
deviation render extremely difficult the determination of the resulting parameter deviations of the elements. 
This difficulty is significantly increased when one takes into account the fact that the reversible and irreversible 
processes of parameter variation may be highly correlated. Under these conditions, an accurate determination 
of the resulting deviations is possible only with an accurate accounting of all the complex disturbing factors andof 
the manifold structural scheme by which the partial parameter deviations are summed when the individual dis- 
turbing factors are operative, 


For an approximate estimate of the resulting element parameter deviations, the solution of the problem can 
be simplified by limiting consideration to the minimum number of those disturbing factors which most significantly 
affect parameter variation, 


In the simplest case, when the sums determining the components of the resulting deviations are made up of 
negligible and mutually independent terms of the same order, all with time-invariant distribution characteristics, 
and there are sufficiently many of such terms, then the distribution law of each of the components may be 
approximately expressed by a normal distribution. 


If the deviations of the disturbing factors are normally distributed then, independently of their number, the 
resulting deviations of the element parameters are also normally distributed. In these cases, the mathematical 
expectation and mean-square deviation of each of the components will be determined by the resulting (total) 
values of the mathematical expectations 


M [Av] = (8) 
M [As] = >) M [yn AAg] 
and total mean-square deviations | 


San = Vt a, {AA,]. 


San = V 28%, o? [AA], (9) 
Sas = V [AA]. 
In this case, the distribution law for the components will be 
{ (B—a,)? 
1 v—a,)® 
v v 
V dno, *P 20°, 


where apg = M [ABy), a, = M [Av], ag = M[ Ae) and Og, Oy and oO, are parameters of the norma! distribution 
functions of the component of the resulting deviations, determined in accordance with (8) and (9). 


In this case, the probability characteristics of the resulting element parameter deviations will be deter- 
mined by the resulting mathematical expectations and the resulting mean-square deviations of these components. 
If the probability characteristics of the components are independent and time-invariant, the mathematical 
expectation of the resulting deviation at each moment of time T may be determined from the expression 


M (AB, (7')} = M [AB,) + {M [Av] + + (M [As] + #5) T?/2, (10) 
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where the values of M [AB ], M [Av] and M [Ae] are 
determined in accordance with (8), 


The resulting mean-square deviation of the ele- 
| ment parameters at each moment of time T will equal 


& 


[AB, (7)] = V cap + (Sav Ty (Sac T?/2)? , 


where the values of Gay and are determined 
sb hc in accordance with (9). 


: Figure 1 illustrates the variation of the resulting 
Time, hours mathematical expectation and mean-square deviation 
of the element parameters under the action of various 
disturbing factors arising in production and use at 

various moments of time (T» through T,). 


No. of defects in 100 hrs, 


Probable Time of Proper Element 


Functioning With Respect to a Given 


Parameter 


If the resulting deviations of the element para- 
meters are random functions of time, then the time of 
proper functioning (useful life) of the element with 

RT ar nr a respect to each parameter can be considered as the 
600 «=6—800— ss inverse functions T = f [AB, (T)] of the given random 
Time, hours argument AB, (T). 


Fig, 6. If the resulting element parameter deviations 
. are given in the form of (6), then the time of proper 
. element functioning with respect to a given parameter is determined by the inequalities (Fig. 2): 


< AB, (1) = ABy + + Aw) T + Ac)T%/2. 


If ABy— AB, arn + (Vo + Av) T + (€ + Ae)T*/2 =0 then, solving this equation with respect to T, we will 
have 


No. of defects in 25 hrs. 


AV (v0 + Ao $2 + As) (AB, — 
e+ Ac (12) 


Expression (12) shows that the time of proper functioning of an element with respect to each of its parameters 
: is a function of many arguments. 


For random values of the arguments, the distribution law of the function g(T) may be determined if the 
distribution law of the arguments and the scheme for forming the resulting parameter deviations are known. 


: In the general case, the computation of the distribution law of a function of several random variables is a 
; sufficiently complicated and difficult task. For an approximate determination of the probable useful life of an 
element it is possible, in this case, to use a simplified computational method. 


If we assume that during the time Tp of proper element functioning the parameter variations may be 
expressed in terms of the mean velocity vme 4nd, further, that with a limited dispersion of the original deviation 
AB, the influence of this dispersion on the resulting parameter deviations can be ignored, then we shall have 


By taking the admissible deviations of the element parameters equal to unity (Badm = 1) and expressing 


the average speed of parameter variation as a fraction of the admissible deviation, we may write Eq. (13) in the 
form 


Tp *® 1/ 


600 
| 
| 
| 
| Umd's < ABadm (13) 
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where Tp is the time of reliable element functioning 
(the reciprocal of the single argument v,,,.). 


If the function v,,¢(T) is continuous and monotonic 
then, based on well-known formulas of probability theory, 
the distribution law of the function ¢(T) of the random 
variable Vie is determined, in the most general form, 
by an equation: 


= %1[%md7)), (14) 


dv me 


where %[Vme(T)] is the distribution function of the argu- 
ment Vines 
>, Since T (Vme) is the reciprocal of V;ne, the distri- 
of T, ¢(T), is determined by the following 


= (15) 


If the argument is distributed normally, (15) takes 
the following form: 


{ 2 


= 


Relative number 
of defects 


Investigation of (16) shows that the distribution curve 
¢(T) which is described by this equation, within the time 
limits - © < T <@, is a symmetric bimodai curve with 
the value zero at T = 0, with two maxima at Ty = 1/¥ 20, 
and T =— 1/ ¥90., sharply ascending towards the maxima 
and then gradually stoping down, approaching the axis 
of abscissas asymptotically (Fig. 3). 


A direct substitution of T = 0 in (16) leads to an 
indeterminate value, but a limiting process shows that 
g(T)~ 0 as T~ 0, since the exponential function 
increases more rapidly than the simple power of T. 


Equation (16) defines the probability density of the 
time of proper functioning of the element with respect to 
a given parameter within the limits ~ «<T < = 


Given the probable useful life of the element, it 
Fig. 9. is of practical interest to determine the probability of 

element failure during time T, counting time from the 
moment of element manufacture (T = 0). 


The probability that the useful life of the element (with respect to a given parameter) lies within the limits 
0 and T or, what amounts to the same thing, the probability of element failure in the time interval T, equals 


P(T)= exp [- dT. (17) 


The magnitude of this integral is graphically defined by the hatched area on Fig. 4. 


8/7) 
4, § 

4 

| 
| 
! 

mit) a 

| 

| 

| 

| 
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The value of the integral in the right member of (17) can, by successive changes of the variable of inte- 
gration, be obtained from the tabulated values of the function of form: 


1 


}=1-—-(z), where z (17a) 


Analogous distribution functions can be obtained for other distribution laws of the disturbing factors and of 
the resulting deviations of the element parameters. 


If the variations of the element parameters are affected, not only by the objective production-usage factors, 
but also by subjective factors (errors of the using personnel), then the probability of element failure with respect 
to each of the parameters will be defined by the expression 


P(T) =1—[1 —Pop (7)| 


— 1, 
or, approximately, 


P(T) = >) Pop (7) + 


where Pop (T) and P,,,(T) are the probabilities, respectively, of element failures due to the objective and sub- 
jective disturbing factors. 


If the probability of element failure due to objective factors is given in the form of (17) but the probability 
of failure due to subjective factors is expressed by an equal-probability law then, with independent disturbances, 
the composition of the distributions will have the form: 


P,(T)=1- (4 ( —cT), (18) 


where c is the parameter of the equal-probability law. 
In this case, the probability of proper element functioning during time T is determined from the equation: 


1 2 
Ry (1) =(1— ( exp [— ar) (1 (19) 


If proper functioning of an element is characterized by n independent parameters where the probability of 
failure with respect to each parameter can be given in the form of (18), then the total probability of element 
failure with respect to all the parameters is, according to (1), 


P,(T) =1— [14 — Pa, (T)I, (20) 
or, approximately, 
Py (7) = >) Pa, (7), 


where the Pa, are determined in accordance with (18), 


The probability of proper element functioning with respect to all its constituent parameters is given in the 
following form: 
(7 


n T 
Re(T)= [] Ro, (7) = n(! \ exp [- 1) (1 


| 


The integrals in the right member can be reduced in accordance with (17a) to the tabulated function cited 
above. 


In the more general case when element malfunction with respect to the individual parameters is char- 
acterized by different probability distributions, the probability of proper element functioning will be determined 
by a more complex composition of the individual distribution laws. 


The analytical expression for the distribution law of element malfunction (16) allows a significant increase 
in precision over the experimental curves obtained by statistical means, 


If element malfunction is expressed in terms of numbers which characterize the ratio of the number of 
failures (m) in a fixed time interval AT to the total number of inspected elements (N), then the quantity g,=m/N 
will approximately express the probability density of improper element functioning in the time interval AT. In 
this case, the curves for malfunction distribution will graphically represent the differential distribution law ¢»(T) 
for malfunctioning. 


As an example, we consider the curve for the distribution of electron tube failures (Fig. 5) which was given 
inref. (3). The given curve was obtained experimentally by tests on the useful life of large numbers of electron 
tubes, a count of failures being made every 100 hours. A consideration of this curve shows that the greatest 
number of failures occurred during the first 100 hours. In each successive time interval, the number of failures 
decreased. The average useful life of a tube turned out to be 583 hours, although the useful life of individual 
tubes reached 3000 hours. 


The distributions of failures of electron tubes differs significantly from a normal distribution, and closely 
coincides with the descending arm of the curve described by relationship (16) (Fig. 4). Registration of the tube 
failures over a comparatively long time interval (100 hours) does not allow a sufficiently accurate representation 
of the ‘failure distribution in the initial time period, from T = 0 to T = 100 hours, to be obtained. If the iube 
failures in the initial period were registered after shorter time intervals (10 to 20 hours), it could be assumed that 
the initial portion of the experimental curve of failure distribution would have the form given in Fig. 6.* 


In this case, the total curve for the distribution of failures in the time interval from T = 0 to T = 3000 hours 
is described, with sufficient accuracy, by Eq. (16).** 


From the physical point of view, the formation of the asymmetric shape of the failure distribution curve 
with the clearly expressed shift of the maximum to the origin is explained as follows: If proper functioning of 
the tubes is characterized by the maintenance of their parameters (tube emission current I) within the established 
admissible limits, then parameter variation, as affected by various usage factors, will be exponential (Fig. 7), 
and these curves may be represented approximately by 


B(T) = B,exp(—a7). 


If there is dispersion in the speed and acceleration of the variation components, the resulting element para~ 
meter deviations in time can be given approximately in the form: 


B(T) = By — (v9 + Av) T + (¢ + Ae) 7/2. 


This function is given graphically in Fig. 8,a. 


The intersections of the curves which characterize the variation of the resulting parameter deviations in 
time with the line parallel to the axis ofabscissas and characterizes the admissible parameter deviation, 
define the limits of dispersion of the time of proper element functioning. 


*Such an assumption is based on the fact that prior to the experiment, the tubes were tested and assorted by the 
OTK with respect to all the parameters, i.e., for T— 0, ¢ (0) * 0. In this case, the number of defective tubes 
in the initial period of the test would be negligibly small, 

* *If, prior to the tests on useful life, the elements had, during a period of time Ty, been stored or subjected to 4 
series of usage disturbances then, as a result of the variations of their parameters during time T,, the elements 
would have had a failure distribution significantly different from the curve described by Eq. (16). 


If the elements had been subjected while stored to the same complex of disturbing factors as those present 
in the test, then the failure distribution laws would have been identical, 


| 
| 
| 
4 
| 
| 


Fig. 10, 


In this case, the distribution of the time of proper functioning of the tubes, i.e., the time during which 
the parameters are maintained within the established admissible limits, is approximately represented by the 
asymmetric curve (Fig. 8,b) which corresponds to the probability distribution law given in (16). 


Dependence of the Law of Element Failure Distribution on the Disturbing Factors 


Given an accurate estimate of the probability of element failure, it is of practical interest to account for 
the dependence of the law of failure distribution on the disturbing factors which characterize the actual con- 
ditions of manufacture and usage of the elements. 


If the probability of element failure is given as the integral function P»(T),then the differential distribution 
law for the probability of failure at the moment T is given in the form 


d|P,(T) 
= (21) 


If the value of the function P,(T) is expressed in terms of the probabilities of element failure with respect 


to each of its parameters [cf. (20)] then, ignoring the influence of the subjective factors, we may rewrite expression 
(21) in the following way: 


Go(T) (1 — (1 — Py — - [1 — Pa 
When we multiply out, and make the replacements dPj/ dt = g;(T) and 1— Pj(T) = R(T), we get 


R, (7) R, (7) R, (T) 


[] 
i.e., 
n Il R; (T) 
= 
(22) 


In the special case when the random variations of the disturbing factors are distributed normally and the 
functions g, (T) and Ry, (T) take the forms given in (16) and (19), the total distribution law for the probability of 


4 (7) 
om: | 
! 
! 
a 
| 
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element failure (with respect to all its parameters) is represented approximately by: 


T 


k=1 0 


After differentiating the right member, we get 


T*V 2no,, 205, 


If, by changing the variable of integration, we bring the integrals in the right members of (23) and (24) to 
the tabulated functions, then (24) will take the following form: 


(T) = ® (z)) 


SUMMARY 
The probability of element failure is determined by the distribution law of its malfunctions, ¢»(T). 


If the functioning of the element is characterized by several independent parameters, the law of mal- 
function distribution is determined in accordance with (22) as the composition of the individual laws of mal- 
function distribution with respect to each of the element parameters ¢j(T). 


The form of the overall distribution law, in this case, is determined by the relative orientations and forms 
of the individual laws, i.e., the functions g;(T), given approximately in the forms of (14) and (16). 


Given an experimental determination of element reliability, the dependence of the overall malfunction 
distribution on the normalized individual distributions in the general form may be expressed by the equation 


malt) 
Nn “| 


where mg (t) = 2 mj (t) is the number of malfunctioning elements, with respect to all the parameters, in a fixed 


n 
time interval At, Nog = > Nj is the number of elements tested, N; is the number of malfunction with respect 
i=-1 
to the parameter Bj, N;/ Ng is the normalizing factor for each individual malfunction distribution and m, (t) is 
the number of malfunctions with respect to the parameter B, in the fixed time interval At. 


Figure 9 illustrates the dependence of the overall malfunction distribution on the normalized individual 
distributions in the case when element functioning is characterized by three independent constituent parameters, 


When the admissible variations of the element parameters are given, the form of the individual distribution 
laws depends on the probability characteristics of the components of the resulting deviations of each of the para- 
meters [cf. (10) and (11). 


The curves of Fig. 10, a illustrate the variation (deformation) of the distribution law of element mal- 
function with unlimited dispersion in the speed of parameter variation and different probability characteristics 
of the resulting parameter deviations, 
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The curves of Fig. 10,b illustrate the variation in the form of the malfunction distribution law with bounded 
dispersion in the speed of parameter variation. 


The probability characteristics of the components of the resulting deviations of each of the element para- 
meters M [4B,), M[v) and M[€],0 ag, and in their turn, depend on the probability characteristics of 
the disturbing factors and are determined in accordance with (8) and (9). If the individual distribution laws of 
element malfunction (with respect to each of its parameters) are obtained,with account taken of the probability 
characteristics of the disturbing factors, then relationship (22) establishes a completely determined dependence 
of the overall distribution law of element malfunction on its parameters and on the probability characteristics 
of the disturbing factors, and allows an analytic calculation to be made of the overall malfunction distribution 
law in various conditions of element manufacture and usage. 


In this case, the limits of admissible deviations of the disturbing factors are determined by the condition 
that a linear dependence be maintained between the variations of the components of the resulting element para- 
meter deviations and the variations of the disturbing factors: 


AB, = Avo = Aey = 


Under these conditions, relationship (22) may be used, not only for the analytic estimation of the pro- 
bability of element failure, but aiso for the theoretical foundations of technological processes and for the regimens 
of accelerated testing of elements for useful life wherein the actions of the disturbing factors are forced. 

LITERATURE CITED 


{1} V. 1. Siforov, "On methods of calculating the reliability of systems which contain large numbers of 
elements,” Izv, AN SSSR, Otd. Tekhn. Nauk, No. 6 (1954) [in Russian}. 


(2) Sh. L, Bebiashvili, "On the theory of the reliability of systems containing large numbers of elements,” 
Izv. AN SSSR, Otd. Tekhn. Nauk, No. 10 (1955) [in Russian]. 


(3) B. Bromberg and R. Hill, Reliability of Aircraft Radio Apparatus Components. Miniaturized Radio 
Apparatus. Translated into Russian. Ed. V. 1. Siforov. Izdatelstvo Inostrannoi Literatury, 1954, 


[4] S. M, Bragin and M, M. Markosian, "The effect of humidity on the electrical characteristics of 
insulating resins " [In Russian), Elektrichestvo, No. 3 (1955). 


(5) V. N. Kalitvianskii, "General regularities in thermal curing of polymerized dielectrics? 
Elektrichestvo, No. 3 (1955). 


[6] V.N. Kalitvianskii and A. V. Kovel’skaia, "Useful life of a new form of electric isolation circuit," 
{in Russian) Elektrichestvo, No. 10 (1955). 


[7] B. Solontsev, "The Use of Piezoquartz Resonators and Exciters [in Russian]. Sviaz'izdat (1955). 


[8] "Some data from the investigation of the basic characteristics of type VS resistors," [in Russian], Inf. 
Sbornik NII MPSS, Nos, 2013-4 (1951). 


(9] F. E. Evteev and V. A. Zhukov, The Technology of Radio Apparatus [in Russian). Gosenergoizdat (1952). 


(10) A. F. Astafev, “Meteorological factors, their effect on radio equipment, and standard methods of 
testing them correspondingly,® [in Russian], Biuro Tekhn. Informatskii MPSS (1948). 


[11] The Reliability of Radio-electronic Equipment (Collected papers) [in Russian]. Sovetskoe Radio (1958), 


[12] I. V. Dunin-Barkovskii and N. V. Smirnov, Theory of Probability and Mathematical Statistics in 
Technology, [in Russian]. Gostekhizdat (1955). 


[13] N. A. Borodachey. Analysis of the Quality of Productive Accuracy [in Russian]. Mashgiz (1946). 


(14) N. A. Borodachev, Basic Questions in the Theory of Production Accuracy [in Russian]. AN SSSR 
(Institute Mashinovedeniia, Leningrad) (1950). 


Received May 30, 1957 


1034 


ON THE NUMBER OF RESERVE SECTIONS 


Vv. G. Druzhinin 


(Moscow) 


Selection of the optimum number of reserve sections in auto- 
matic systems is examined, 


The use of reserve (spare)elements is lately becoming of increasing importance, This is caused by auto- 
mation becoming more and more complex and by the increase in the number of elements used; this has not been 
matched by a corresponding improvement in their quality. A number of papers [1-4]has been devoted to the 
investigation of various reservation problems. Of greater practical importance is the problem of reserving in 
systems which have switches automatically cutting off the defective elements and switching in the reserve ones 
[2, 3]. It has been shown in our literature that the reserving of separate sections (called sectional reserving) is in 
most cases more convenient than the reserving of a system as a whole (called overall reserving) and also shown what 
should be the lower bound of the reliability of a switch so that sectional and overall reserving are equally effective. 
The problem, seen in this way,is certainly ofinterest. From practical experience, however, a number of other 
problems have arisen which so far have only been investigated a little. One of these is briefly considered below. 


2 


OF 


Fig. 1, Reservation scheme. A) basic and reserve 
elements, B) switch circuits, m) number of reserve 
sections, k) number of reserve elements per 
section including basic elements. 


Fig. 2. Probability Q of breakdown of a reserve 
system with different numbers m of sections when 
k = 2 and qs = 0.01. The dependence of Q(m) on 


The engineer is usually confronted with various 
limitations when constructing a system with reserve 


elements and with switches (afterwards called reserve m for various Q, is calculated from formula (4). 
systems). The multiplicity of reserving, that is,the Crosses denote points of the lowest value of function 
specified number of sections which have reserve units, Q(m). 


is usually limited by reason of cost, weight and volume. 

The reliability of the switches depends on the reliability 

of their components, on some characteristics of the system performance, and on other factual points. Therefore, 
the problem the engineer is faced with is,in the first instance, as follows: into what number m of sections should 
one divide the whole basic system which has probability of breakdown equal to Qg™ so as to obtain the highest 
possible reliability, which means the least possible probability Q of breakdown of the system? We denote the 
probability of breakdown of a switch by q, and the total number of units in each section by k. When the problem 
is thus formulated, the overall reserving appears as a particular case of sectional reserving with m = 1. 
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Fig. 3. The least probability of breakdown of the reserve system as dependent 
on Q, for several q, when k = 2, 
We shall assume that the main and reserve circuits have 
the same performance capability. 


The probability qp of breakdown of one section of the 
basic network is connected with the probability Q, of breakdown 
of the overall basic network as follows: 


Qo = 1—(1 — (1) 
Whence 


1 
qo=1—(1—Q)™. (2) 


| Using (2) one can write the formula for the probability of breakdown of one knot, the knot being one 
section with all its reserve units, as: 


11 — (1 —9)(1 — 


As the overall reserve system (Fig. 1) consists of a succession of knots with reserve units, the probability 
Q of its breakdown is given by the expression 


1 
Q=1—({1— [1 —(1 (1 — ]*)™. (4) 


We shall consider systems in which distribution functions Qy (t) and qs (t) take only small values during a 
comparatively short operation time of the system. If so there exists such a number of sections that the probability 
of breakdown of the overall system becomes least (Fig, 2). The most suitable number of sections m will be 
found when Q, « 1 and q; « 1 for a specified operation time T. 


The right-hand side of (4) is successively expanded in powers of 1/m and m; then the expressions in paren- 
theses are multiplied out,neglecting the terms of low order of magnitude, and hence we obtain 


Q =m [2s +- Qo} . (5) 


4 
| 


The maximum and minimum values of function Q (m) as given in (5) are now found, The solution ofthe 
equation 


aQ k 
gives the result m =mythe extremal value of the probability Q of the breakdown of the system. This value fs 
my (kA — 1). (7) 
Substituting the value m = my into the second derivative a*Q/ i) m*, one can verify that the function Q (m) 


as given by (5) has at the point my a minimum. 


Thus for small values qs and Qo, the least probability of breakdown of a reserve system, and therefore the 
highest efficiency of reserving, occurs when the number of sections is equal to mysgiven by formula (7). The 
values of my should be rounded off to the nearest integer, 


The accuracy of formula (7) may be estimated from the table in which the values of the percentage error 
are only given, for k = 2, this being the most important case from practical point of view. 


The percentage error is evaluated from the formula 


— my, 
100%, 


m ex Mex 


where m,, is the exact value of the optimum number of sections obtained from the exact formula (4) for Q (m). 


The increase in error for higher values of the ratio Qo/qs is not of great importance as in this case the graph 
of Q (m) is very flat in the vicinity of the minimum value (Fig. 2). 


Substituting into (5) the value of my frum (7) we oltain an expression for the least probability of breakdown 
of the system 


Ke 
Qmin= (8) 


Formula (8) allows estimation ofthe reserving efficiency when Qo, qs, k are given. 


Figure 3 shows the dependence of the least probability Q,,jn, Of breakdown of the system on Qy for several 
possible values of q,, when k = 2, The continuous line denotes the relationship Qrin (Qo) obtained from the curve 
Q (Qy ds, m, k) according to formula (4) the dashed line represents Q,,j,, found from the approximate formula 
(8). 


Expanding the right-hand side of formula (2) by binomial expansion and neglecting powers of Q, higher 
than first, the probability of breakdown of one reserve section can be written in the form 


qo = Qo (9) 


m 


Substituting into formula (9) the value of my from (71) we obtain the relationship 


(asopt = (k= (10) 
If follows from (10) that reserving attains its highest efficiency when probability q, of breakdown of a 


switch is proportional to probability q9 of breakdown of one section, and the coefficient of proportionality appears 
as the number of reserve units (k~ 1) in each section, 


The formulas obtained above are true for fixed periods of time; but if the condition Qg (t)/q, (t) = const is 
fulfilled, then they are true for any period of time during which the system is operating. 


If Qo (t) and q, (t) are two different time functions,then the optimum number of reserve sections as evaluated 
from formula (7) may somewhat vary with the length of time of operation. In this case an average value of my 
over equal periods of operation of the system is found, 


Am _ Mex 
| 
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Qo 
0.01 | 0.05 | 0.40 0.20 | 0.30 | 0.40 
0.01 | O 0 0 0 |4-9% [420% 
0.10 0 0 0 0 0 0 


When choosing an average value of m, one should take into account that for higher values of ratio Qo/4q, 
the error in my is not importantgs was shown above. 


Formulas (7), (8) and (10) may be of great help in finding a practical solution to this problem. One should, 
however, take into consideration that,in practice,the number of reserve sections is determined not only by their 
reliability, but also by a number of design factors which often prevent the full application of the obtained formulas. 
Also, in actual use, the reliability of reserve sections and of switches may vary from one section to another. There- 
fore, the results obtained from the above formulas should be regarded as tentative. 
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EFFICIENCY OF USING FREQUENCY BANDS IN TELEMETERING 


R. R. Vasilev 


(Moscow) 


A method is considered of determining the efficiency of using frequency bands 
in telemetering channels for various types of modulation. The method relies on the 
application of information criteria in the evaluation of telemetering systems. 


A telemetering system is used to transmit information on processes which occur at the transmitter. To 
evaluate the performance quality of such a system it is necessary to determine what quantity of information on 
the measured processes can be obtained at the receiving end of the system. The telemetering device including 
transmitters may be considered as an information source, and knowing the parameters which determine the per- 
formance of the system,one can evaluate the rate with which the information will be generated at the receiving 
end of the system, Other conditions being equal, the system which provides more information is more precise. 
The quantity of information can be measured, and its magnitude provides the basis for evaluation criteria of the 
performance of a telemetric system. 


A number of criteria for evaluation of telemetering devices was given in paper [1)}, and three different 
systems in actual use were also compared, In the same way, other systems not yet in use can be investigated, By 
applying information criteria one can estimate the advantages of various systems of signal modulation in tele- 
metering devices. To simplify calculations we shall consider single-channel telemetering devices only, with 
different types of signal modulation; but the method in question is also suitable for multi-channel systems, The 
system evaluation will be carried out from the point of view of efficiency of using frequency bands in communi- 
cation channels, Telemetric systems can also be investigated, by using analogous methods, from another point 
of view (for example cost, complexity, reliability etc). 


» Ry, bits 
AM 2/¢ (logs ep — 2.5) / ¥ 
FM 0.16 / p? (y — 1)* [logs e (y — 1) — 0.68] / 
PFM 0,1(1—3 / y + 9.2 /y*)/ 
PTM 0.0625 / (logapy'*) / 
PoM-t | 0,125 / (logepy'* — 0.5) 
PCM-2 0,25 / (logspy'* — 1) 
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Fig. 1. Dependence of Rg on p for AM (continuous 
line) and FM (dashed line). 
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Fig. 3. Dependence of R¢ on p for PTM (continuous line), 
for PCM-~-1 (dashed) and PCM-2 (dash and dot). 
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Fig. 2. Dependence of Rr on p for PFM. 


It is assumed that the measured quantity 
has a rectangular probability density distribution 
in the interval (— a,a). Spectrum of the measured 
quantity is in the frequency band 0-F,,. Thus, 
the transmission of different values can be carried 
in time interval T = 1/2 F,,. The measured 
quantity is recorded at the receiving end,to- 
gether with error due to noise in the communi- 
cations lines. One can assume the error to be 
distributed normally (which is correct for most 
types of modulation when noise is weak), To 
evaluate the error in such cases it is convenient 
to introduce the magnitude of relative error 


3’ 


where §* is the mean-square error and 2a the 
variability range of the measured quantity. 


The output rate R of the messages at the 


exit end of the telemetering system is, in accordance with the assumptions made, bounded [1] by 


2 
Fm logy < R< Frm logs (2) 


The inequality can be written by using the introduced relative mean-square error. 


From (1) and (2) we obtain 


Fm loge Fm logs (3) 


The upper and lower bounds of R differ by approximately 0.5 F,,. Therefore, with an accuracy sufficient 


for practical purposes 


1 


Further on we shall replace the approximation sign by ordinary equality sign. To evaluate the efficiency 
of using frequency bands in communication channels it is essential to apply the R¢ criterion which gives the quantity 
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4 
R = Fm logs (4) 


of information per unit frequency of communication channel. This criterion is given by the formula 


R_ | 


where Af is the frequency band of the communication channel, 


The magnitude of 6 can be determined by using methods of potential interference rejection. In paper [2] 
values of 5 are given for different types of modulation with weak noise and with peak-limiting of signal in the 
communication lines. In the paper the magnitude 6 is expressed in terms of generalized parameters 


U. 
= 
where U, is the signal amplitude in the communication line, 09 the noise voltage. 
The table gives values of 6" and R, obtained from formula (5). 


Figures 1-3 show graphically the dependence R¢ on p, that is Rg = f(p). 


When the signal is amplitude-modulated, it is sufficient to have the channel band Af equal to 2F,,, (y = 1). 
In many cases, however, when Fry is small and only limited facilities available,one is unable to obtain a suf- 
ficiently narrow band of communication channel, that is to have equality Af = 2Fm- For this reason function 
Ry = f(p) is shown in Fig. 1 for various y. Figurt 2 shows that in the case of PFM, for different noise levels in 
communication channe! there exists an optimai value of the coefficient of band broadening y = Yopt. The 
analytic determination of y op, becomes difficult, as the equations leading to it are cumbersome, 


Figures 1-3 show that in the region of very weak noise the amplitude modulation when y = 1 secures the 
best use of the frequency band of the communication channel, If, however, one is unable to obtain a sufficiently 
narrow communication channel band necessary for AM, and one has to operate in the region y» 1, then, in the 
given meaning,the amplitude modulation becomes inferior to other types of modulation, 


FM, PFM, PTM and PCM have approximately equal efficiency in frequency band usage, especially for 
greater values of y. 


If in the determination of 6 one takes into account not only noise in the communication line but also 
instability of equipment, one is able to acquire a better appreciation of the efficiency of using communication- 
channel frequency bands. Unfortunately, no detailed data on instability of equipment is available in the literature, 
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STABILIZATION OF TEMPERATURE OF OPERATING SUBSTANCE IN THERMISTORS 


P. Udalov 


(Moscow) 


A diagram is given of temperature stabilization of operating substances in ther- 
mistors with accompanying temperature changes in the surrounding medium, as well 


as the calculations for the setup. 


The use of semiconductor thermistors (STR) as inertial elements in time relays, as well as in many other 
electrical devices which are based on nonlinearity of the volt-ampere response of STR, is limited by high sensitivity 
of such devices to oscillations in temperature of the ambient medium, However, by using an indirectly heating 
STR one is able to construct a simple device (Fig. 1) which ensures temperature stabilization of the STR operating 
substance by means of changes in the heating current that will depend on the ambient temperature. 


Fig. 3. Family of volt-ampere characteristic curves of 
bypass STR. 


4 


Fig. 2. Heater characteristic of stabilized STR 


The device operates in principle as follows. 
With the decrease of temperature 6 of the medium, 
current }}, in the coil of the heater increases as the 
magnitude of resistance Rpp of the bypass STR goes up, 
but current I, in the linear circuit of the network is 
kept constant by ballast resistor Rp. 


Instead of Ry one can also utilize a linear 
resistance whose magnitude is so chosen that the change 
in current Ip is negligible as compared with the change 
in current I}, in the heater coil. In this case, however, 
a considerably higher voltage is required than in the 
case of a ballast resistor. 


In order to carry out the calculations for the 
device in question, it is necessary to have the heating 


characteristic curve of the stabilizéd STR [1], which expresses the heating effect of the heater current in the STR 
circuit on the operating substance (Fig. 2); it is also necessary to have the family of volt-ampere characteristic 


curves of the bypass thermistor Rpp (Fig. 3). 
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Fig. 1. Stabilization network of temperature of STR 
working substance. 
% min 
6,, 
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A, 


‘h min . 


Fig. 4. Determination of Rh, Ig and Ih = f(@9): a) 


volt-ampere characteristic curves of bypass STR, 
b) heater characteristic curves of stabilized STR. 


Fig. 5. An example of graphical construction. 
When the temperature @, of the medium varies 
from 99min tO 89max then the temperature @7 of the operating substance should remain constant and equal to 


= + Oh = Somax- 


This relationship can be made true by yarying the second term from @ (when @» = @gmax) to A@g = Ogmax ~ 
— @omin- This can be realized if the current in the heater changes from zero to i, depending on the temperature 
of the ambient medium in accordance with the graph of heater characteristic. But at Ogmax the resistance 
does not become zero, and so the zero value of heater current I, and of temperature 6} cannot be reached, 
Therefore, at temperature @gmax there will be some current Ip min in the heater, and the temperature 6-7 of the 
operating substance will be higher than O94, The value of Ip min depends on the actual form of heater char- 
acteristic and on the selection of working section AB (Fig. 2). 


The selection of a working section of heater characteristic is determined by two factors: the section should 
be as nearly linear as possible, and temperature @} corresponding to current Ihmin should be minimum so as not 
to cause redundant overheating of the operating substance. 


The following graphical method can be used both to determine resistance Rp of the heater coil and current 
I, in the linear circuit, as well as to obtain the characteristic , = f (@9) which, as a rule, is not identical with 
the selected section of heater characteristic. 


Volt-ampere characteristics of the bypass STR Rpp and the heater response of the stabilized STR Ry are so 
drawn that the lines representing currents are parallel (Fig. 4). From end points A and B of the working section 
of heater characteristic two lines are drawn parallel to the ordinate axis; their points of the intersection A’ and 
B’ with volt-ampere characteristics should in our case be collinear with point C which is the point of inter- 
section of the ordinate axis of heater with the abscissa axis of volt-ampere characteristic (Fig. 4). 


If it is impossible to find one straight line through points A’, B’ and C, a straight line is drawn through C 
equidistant from points A' and B*, The abscissa of point C is the required value of current Ip, and the straight 
line itself is the volt-ampere characteristic of resistor Ry, of linear circuit whose slope is a = tan~* Rp. 


* 
! 
| 4 
| 
“hmax 
h 
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The points of characteristic curve I, = f (69) can be obtained by projecting the points of intersection of 
the linear circuit volt-ampere characteristic curve with the Rbp volt-ampere characteristic curve onto straight 
lines drawn parallel to the current axes through the corresponding points of the ordinate axis of the heater char- 
acteristic curve, The characteristic curve thus obtained shows the actual dependence of the heater current on the tempera- 
ture of the medium in the section under consideration.It may occur, however, that the curve I) =f (@ 9) thus obtained 
will be wholly above or below the working section of the heater characteristic curve, In this case, by varying the slope 
angle a of the volt-ampere characteristic as well as the magnitude of the current Ip previously found, then by trial 
and error the curve Ih = f (@9) can be placed in such a position with respect to the working section of heater 
characteristic curve, that the error becomes least. 


The task becomes much easier if one of the graphs, say heater characteristic, is plotted on tracing paper, 
and the other on graph paper. 


To determine the voltage of the current feed U at point C, the volt-ampere characteristic of the ballast 
resistor is drawn, as shown in Fig. 4. When instead of Rj, a linear resistor is used, the current Iy will not remain 
constant. To reduce the error thus introduced ,the resistance Rp is so chosen that the computed value of current 
I, occurs when the voltage drop in the heater (Fig, 5) is equal to 


Uhmin + Ub max 
Up = 3 


The proposed design has been verified experimentally, for temperature @9 varying from 0 to + 50°C. 


Temperature stabilization was attained of STR working substance, type TKP-450; its heater curve 
is shown in Fig. 5,b, and the temperature response is expressed by the formula 


4750 
Ry = 23.7 x 10%, T ohm. 


Two STR, type TOS-M, connected in parallel, were used as’ Rp. The family of their volt-ampere character- 
istic curves is snown in Fig. 5,a, and the temperature response is expressed by the formula 


2900 
Ry =0.143e T ohm. 
Resistance Rp is linear. Power voltage U = 220 volts. 


The operating section of heater characteristic was so selected that 6}, was between 20 and 70°C; thus the 
temperature of the operating substance 6-7 should have been kept constant by the stabilizer at 70°C. 


Resistances Ry, = 10620 ohm (Fig. 5) and R,, = 918 ohm were obtained. 
The results of measurements of resistance in the stabilized working substance, and the corresponding 
temperatures 6-7 are given in the table below: 


°C 0 
Rr, ohm 2100 
Or, 714 
SUMMARY 


1, By introducing the outlined network for the stabilization of temperature of the operating substance, the 
performance of time relays, generators and of other apparatus with STR can be improved considerably. 


2, The graphical method of computation of the stabilizer’s circuit diagram gives satisfactory accuracy for 
practical purposes, 
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LOGARITHMIC METHOD OF PLOTTING REAL FREQUENCY RESPONSE 
OF AN AUTOMATIC CONTROL SYSTEM 


Vv. A. Atsukovskii 
(Moscow) 


A logarithmic method of plotting real frequency response of a system is proposed 
which does not involve cumbersome calculations, The method also permits one to find real 
roots of algebraic equations of any degree. 


The transient performance of a closed-loop servo system is often obtained from its real frequency response, 
The disadvantage of this method lies in rather complicated calculations necessary to graph P (w) = ReW (jw). 


The proposed method permits evaluation of roots to three significant figures if plotting is accurate, It uses 
tables of logarithms of sums and differences computed by Gauss, The sum or difference of two numbers can then 
be found from their logarithms. Tables computed by Gauss giving the logarithms of sums and differences are 
available. 


Finding the logarithm of a sum. Given log a and log b, where a > b, a and b of equal sign. The difference 
R = log a— log b is worked out. The 2 corresponding to R is found from the table, and hence log (a + b) = loga +2 


is obtained, where 


cay]: (1) 


TABLE 1 


SOS 9 


ooo 


TABLE 2 


Expressions 
hing pw) values at 


Igo = 0 Igo = 1 


Constant factor 
k = 6,5« 108 


Numerator terms 
1232 
3.69 x 10° 


Denominator terms 
8iga 
6.54 x 10a 3.8164 6lg 
1.3 x 6.112+44]lg 
2.9 «10%? 9.462+42I¢ w 
2.22 10? 12,345 


i 
garithm of root 


Fig. 1. 
The sum will have the sign of its terms. 
Finding the logarithm of a difference, Given log a and log b , a > b, and a and b having different signs. 


The difference R = log a— log b is worked out. The A corresponding to R is found from the table and thus 
log (a— b) = log a— A is found where 


A=lg-———. 
| 

The difference takes the sign of the greater figure. 

A few values of = (R) and A (R) are given in Table 1. 


The values of © and A are tabulated more extensively in [1]. 


Logarithms of 
4lga 0.000 4.000 
2.09-+-2Igeo 2.090 | 4.090 
5.566 5.566 5.566 
0.000 8.000 
3.816 9.816 
6.112 10.112 
9.462 11.462 
12.345 | 12.345 
13 
12 
g 
ae 
7 
6+ + 10 
§ 
4 
2 
0 2 
Lo 
(2) 
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Fig, 4. 


If a sum of three numbers is required then we 
first find the sum of two numbers by looking it up in the 
table; then the obtained sum (or difference) is added 
to the third. 


The general form of the rea! frequency 
response is as follows: 


+ +...4+A, , Mo) 
Bow” + +4. . By 


p(o)=k 


where k is a real,constant factor, A‘s and B's are constant coefficients, a and 6 are degrees of the polynomials 
in the numerator and denominator respectively. 


The logarithms of each term of the numerator and denominator are taken separately; 


Ig Ig Ao + alge, 
lg = Ig Ap + (a — 2) Iga, 


Ig A, == Ig A,; 


ign 
26 
2 
22 3 cous 
2 
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Ig Bow” = lg By +- B lg w, 
Ig Bywo** — Ig Bs + (B — 2) lg a, 


The logarithm of each term becomes a straight line in logarithmic scale, and can be plotted through two 
points log w = 0 and log w = 1. The obtained expressions are plotted separately for the numerator and also 
for the denominator, the values are then added at various points by using the table. Having added the logarithm 
of constant factor to logarithm of numerator and having subtracted logarithm of denominator, we obtain the 
logarithm of the system's real frequency response as a function of w. 


Example, We shall plot system's real frequency response as given by 


6.5 x 10°(—w*— 1230 + 3.69 x 10) 
P + 6.54 — 1.3 x — 2,9 X10? 2.22 108° 


We find logarithms of each term separately, both in the numerator and denominator, The calculations are 
shown in Table 2, 


Figure 1 shows how the graph of the logarithm of numerator is plotted, Fig. 2 of denominator, and Fig. 3 
how the logarithm of the real frequency response is obtained by adding the logarithms of constant factor and of the 
numerator, and subtracting the logarithm of denominator. Figure 4 shows the real frequency response itself 
obtained from its logarithm. 


When plotting logarithms of these polynomials we determine the sign, for a given value of w, by the rule 
given above, and the sign of the whole real frequency response is obtained from the sign of numerator, denominator 
and multiplier for any chosen value of w. 


The method can be used to evaluate real roots of polynomia! expressions. When plotting logarithm of the 
polynomial at point w representing a root, logarithm of the polynomial tends to and the polynomial expression 
itself changes its sign. 
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ANNOUNCEMENT 
The Institute of Automation and Telemechanics of USSR Academy of Sciences announces that the 3rd All 
Soviet Conference-Seminar on Pneumo~Hydraulic Automation will be held in April-May 1959. Reports on the 
first two conferences were published in the journal "Automation and Telemechanics," No. 12, 1957 and No. 10, 
1958. 


The Institute of Automation and Telemechanics of USSR Acad, of Sci. invites research and development 
institutes, drawing offices, factories and specialists in planning, production, maintenance or utilization of pneumo- 
hydraulic automation, to attend the 3rd Conference-Seminar, 


Applications together with lists of participants, as well as contributions or communications with brief 
resumes should be sent to: Moscow, Kalanchevskaia Str., No. 15a, Inst. of Autom. and Telem, of USSR Acad, of 
Sci,, Laboratory of Pneumo-Hydraulic Automation, not later than February 1, 1959, 


Organizing Committee of Conference-Seminar on Pneumo-Hydraulic Automation. 


